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LIGHT  SCAT  TLR1NG  KROM  IRREGULAR  DILLKTRIC  -’AKTIC.l  IS 

Gray  Ward 

Introdjii  Mon 

Accurate  exper i mental  measurements  of  scattering  from  irregular 
dielectric  particles  of  known  size  and  shape  are  needed  to  verify 
theoretical  calculations  of  scattering  from  these  objects  when  the 
wavelength  of  the  incident  radiation  is  comparable  to  the  particle 
dimension.  Holland  and  Gagne  [1|  have  made  careful  measurements  from 
irregular  silica  particles  from  18°  to  166",  and  find  significant 
deviations  from  spherical  models.  Parmenter  L 4 has  measured  the 
angular  distribution  of  scattering  from  randomly  oriented  single 
particles  of  quasi-el  1 ipsoidal  BaSO^,  and  also  finds  significant 
deviations  from  spherical  particle  scattering  near  backscatter. 

Pinnick,  Carroll  and  Hofmann  [5]  measured  the  scattering  from 
clumps  of  NaCl  and  M>0,  particles,  and  found  that  for  particles  with 
size  parameters  larger  than  about  b,  Mie  theory  overestimates  the 
scattering  at  backscatter  angles. 

Objectives 

L I OAR  measurements  of  atmospheric  aerosols  utilize  the  spherical 

particle  model  and  look  at  backscatter.  We  wished  to  extend  the 

measurements  of  scattering  from  irregular  particles  to  angles  very  close 

to  backscatter  to  see  how  serious  an  error  is  being  introduced  bv 

wanted  to  determine  whether  an 
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ensemble  of  tumbled  irregular  particles  could  be  represented  bv  an 
equivalent  sphere  in  the  backscatter  ranqe  as  it  can  be  in  forward 
scatter. 


Approach 

Mie  scattering  from  low-loss  spheres  of  one  size  has  very 
pronounced  structure,  which  is  washed  out  by  size  dispersion  or 
by  departure  from  spherical  shape.  We  elected  to  work  with 
optical  scattering  in  the  laboratory  with  aerosols  of  controlled 
size,  shape,  composition,  and  size  dispersion.  Monodispersions 
of  particles  of  known  size  and  shape  were  generated  with  a 
Llerglund-Liu  [6]  v ibrati ng-ori fice  aerosol  generator.  Particle 
shape  and  size  distribution  was  monitored  bv  SEM  photographs. 

The  light  source  was  a pulsed  nitrogen  laser  pumping  a 
dye  laser.  Pulses  were  typically  10  kW  4 ns  FWHM  at  627  nm, 

10  pps. 

The  polar  nephalomet.er  went  through  three  models.  The  first 
model  used  a scanning  mirror  mounted  on  an  arm  to  sample  and 
reflect  a V'  x 1"  portion  of  the  scattered  light  to  a small  mirror 
on  the  axis  of  rotation,  and  out  through  the  tube  on  which  the 
arm  pivoted.  The  scanniny  was  controlled  bv  a stepping  motor. 

Light  was  analyzed  bv  a fixed  quartz  analysis  prism,  detected 
by  photomultipliers.  Mechanical  overshoot  was  a problem.  The 
second  model,  used  to  obtain  the  data  reported  herein,  used  96 


fixed  mirrors  in  place  of  the  scanning  arm  (figure  1).  Scanning 
was  accomplished  by  rotating  a small  1.2  cm  diameter  mirror  on 
the  axis  of  rotation;  though  controlled  by  a stepping  motor, 
the  system  had  no  problems  with  mechanical  overshoot.  The 
scanning  was  accomplished  in  23  sec,  limited  by  puner  tape  punch 
rate.  The  third  model,  which  is  being  tested  at  present , 
eliminates  scanning  altogether.  It  uses  fiber  optics  at  angles 
selected  by  a Backus-Gi lbert[7]  routine,  of  various  lengths  so 
as  to  delay  each  scattered  pulse  by  a different  amount,  forming 
an  optical  paral lel-to-serial  pulse  converter.  The  longer  fibers 
not  only  have  more  delay,  they  have  higher  loss,  and  they  are  placed 
strategically  at  angles  where  scattering  is  stronger,  so  there 
is  compression  of  the  dynamic  range  of  the  detector  output.  The 
fiber  bundle  goes  to  a photomultiplier  detector  (figure  ?.) . 

The  chamber  in  the  third  model  contains  provisions  for  suspension 
of  a single  oriented  particle,  and  is  vacuum-tight.  We  are 
attempting  to  capture  single  particles  after  scattering  on  a glass 
slide  and  to  prepare  them  for  SEM  photography  in  the  same  apparatus. 
The  third  model  provides  a pulse  train  with  the  scattering  information 
obtained  from  a single  laser  pulse  incident  on  a single  oriented 
particle. 

The  data  acquisition  system  was  designed  and  built  at  the 
University  of  Florida,  which  was  a mistake.  It  delayed  the  program 
at  least  two  years.  It  consists  of  a fast  sample-and-hold  unit  for 
each  of  two  photomultipliers,  one  for  each  polarization  channel,  an 
A/D  converter,  paper  tape  punch,  and  control  logic.  It  could  capture 

the  4 ns  pulse  and  perform  a 12-bit  digitization  reliably.  (The 
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development  problems  were  due  to  noise  pickun  from  the  thyratron  in  the 
nitrogen  laser,  combined  with  personnel  changes  in  the  Electronics 
Shop.)  At  the  time  it  was  built,  it  exceeded  the  capabilities  of  any 
commercially  available  equipment.  For  the  third  model  nepha Icimeter , a 
different  data  acquisition  svstem  is  required,  such  as  Tek  R7912/w»’  <010. 

The  procedure  was  to  build  the  apparatus,  test  with  spheres  and 
compare  with  Mie  theory,  test  with  ellipsoids  and  compare  with 
published  results  of  others,  then  proceed  to  cubes  and  other  shapes 
we  manufacture. 


Exper ime n tal  Resu 1 t s 

The  calibration  of  the  apparatus  is  shown  in  Figure  3.  The  upper 
curve  shows  the  experimental  data  points  with  individual  error  bars 
together  with  the  theoretical  curve  calculated  from  Mie  theory, 
using  the  index  of  refraction  and  size  dispersion  data  provided  by 
the  manufacturer.  There  are  only  two  adjustable  constants  altogether: 
the  gain  of  each  photomultiplier.  Calibration  runs  are  taken  before 
and  after  each  experimental  run  with  non-spherical  particles.  In 
this  case  we  would  expect  the  data  to  agree  with  the  theoretical  curve, 
and  to  the  extent  that  it  does  not,  further  improvement  in  our 
experimental  procedures  is  indicated.  A SEM  photograph  of  the  calibration 
particles  is  included.  The  fit  in  Channel  2 is  very  good,  and  in  Channel 
1 it  is  fair. 
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The  data  taken  from  Cl adospori um  particles  is  shown  in  Figure  4. 

The  theoretical  curves  are  the  best  fit  that  could  he  obtained  by 
using  a spherical  particle  model  and  varying  the  size  parameter. 

The  fit  is  qualitatively  correct  in  the  middle  range  of  scattering 
angles,  but  does  not  fit  the  data  as  well  as  in  the  calibration  runs. 

This  difference  in  goodness  of  fit  indicates  that  the  tumbled 
Cl adospori um  particles  are  only  crudely  represented  by  a spherical 
scatterer.  An  SEM  photograph  of  a cladosporium  spore  is  shown. 

(The  SEM  preparation  process  has  dried  out  the  spore  so  that  it 
has  partially  collapsed,  but  its  original  spheroidal  shape  is 
discernable. ) 

The  data  from  NaCl  crystals  is  shown  in  Figure  5. Again  the  theoretical 
curves  are  the  best  fit  that  could  be  obtained  using  a spherical 
particle  model  and  varying  the  size  parameter  (the  index  of  r 'fraction 
is  known).  The  SLM  photograph  shows  that  we  have  avoided 
clumping  of  the  salt,  crystals  through  beta-particle  irradiation  of 
the  aerosol.  ( It  is  questionable  if  the  particles  are  rounded  on 
only  one  side  as  shown;  the  shape  has  been  influenced  by  the  collect!  'i 
and/or  sample  preparation  process.  The  volume  is  correct.)  Wo  cce 
that  the  spherical  particle  model  is  qualitatively  correct,  but  it 

t 

does  not  fit  the  data  as  well  as  it  did  in  the  calibration  runs. 
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SCATTERING  ANGLE  • . 

Figure  5.  Quasi -cubical  NaCl  Crystals  (4 pm)  with  best  fit  SDherical  model (k=l 0.2) 


Cone! usions 


We  have  constructed  a digital  polar  nephalometer  which  takes 
accurate  detailed  scattering  patterns  from  irregular  dielectric 
particles  of  known  and  irregular  shape  in  tumbled  monodispersions. 
The  accuracy  of  the  apparatus  at  this  point  is  indicated  by  the 
calibration  runs  in  Figure  3.  We  have  taken  subsequent  scatter  inn 
patterns  from  two  different  types  of  irregular  particles  about 
which  we  have  index  of  refraction  information  and  some  shape 
information.  The  data  is  not  well  fit  by  any  sort  of  equivalent 
sphere  model;  therefore  tumbled  irregular  particles  do  not  scatter 
like  eqivalent  spheres. 

Further  work  is  warranted.  The  data  is  of  sufficiently  high 
quality  to  lend  itself  to  fitting  with  theories  for  irregular 
particles[8]-[12]. 


[1]  A. C. Holland  and  G.  Gagne,  Applied  Optics  9,  1113(1970) 

[2]  T.R. Marshall,  C. S. Parmenter , M.Seaver,  Science  190,375(1975) 


[3]  T.R. Marshall,  C.S. Parmenter,  M.Seaver.J.Col loid  and  Interface 

Science  55,  624(1976) 

[4]  C.S. Parmenter,  private  communication. 

[5]  R.G.Pinnick,  D.E. Carroll  and  D.J .Hofmann , Apnlied  Optics  15, 

384(1976) 

[6]  R.N.Berglund  and  B.Y.  Liu,  Env.  Sci.  ft  Tech.  7,  147(1973) 

[7]  E. R.Westwater  and  A. Cohen,  Applied  Optics  12,  1340(1973) 

[8]  G.Ward,  Proc.  IEEE,  Southeastcon  74,160(1974) 

[9]  C.Yeh,  P. Barber,  UCLA-ENG- 7431  (1974) 

[10]  Reiner  Eiden,  Appl ied  Optics  5,  569(1966) 

[11]  Herschel  Weil,  C.M.Chu.  Applied  Optics  1 5 , 1832(1976) 

[12]  S.Asano,  G.  Yamamoto,  Applied  Optics  14,  29(1975) 


12 


SCATTERING  OF  LIGHT  FROM  NONSPHERICAL 
DIELECTRIC  PARTICLES: 
SPHEROIDAL  CLADOSPORIUMS  (2pm) 
CUBICAL  NaCl  CRYSTALS  (4pm) 


ISMAIL  KIRMACI 


A THESIS  PRESENTED  TO  THE  GRADUATE  COUNCIL  OF 
THE  UNIVERSITY  OF  FLORIDA 
IN  PARTIAL  FULFILLMENT  OF  THE  REQUIREMENTS  FOR  THE 
DEGREE  OF  MASTER  OF  ENGINEERING 


UNIVERSITY  OF  FLORIDA 


1978 


TABLE  OF  CONTENTS 


I 


Page 

ACKNOWLEDGMENTS  iii 

ABSTRACT  vi 

CHAPTER  I:  INTRODUCTION  1 

CHAPTER  II:  THEORY  8 

Theoretical  Background  8 

CHAPTER  III:  EXPERIMENTAL  APPARATUS  12 

Scattering  Chamber  12 

Digital  Data  Recording  System  14 

Operation  15 

Nitrogen  Laser  16 

Aerosol  Generator  17 

CHAPTER  IV:  EXPERIMENTAL  RESULTS  AND  CONCLUSION  29 

Corrections  for  Experimental  Data  33 

Experimental  Data  And  Graphs  36 

Analysis  of  Data  And  Error  Calculations  85 

Conclusion  91 

APPENDIX  I:  Electromagnetic  Scattering  95 

APPENDIX  IT:  General  Theory  Of  Scattering  by  A 

Sphere  103 

iv 


Page 


APPENDIX  III:  Scattering  By  An  ellipsoid 


APPENDIX  IV: 


APPENDIX  V: 


APPENDIX  VI: 


APPENDIX  VII: 


Computer  Program  To  Decode  Paper 

Tape  And  Average  Over  Statistical  Data  124 

Theoretical  Mie  Program  132 

Volumetric  Correction  136 

Scattering  Expectations  141 


APPENDIX  VIII:  Data  Tabulation 


LIST  OF  REFERENCES 


BIOGRAPHICAL  SKETCH 


Abstract  of  Thesis  Presented  to  the  Graduate  Council 
of  the  University  of  Florida 
in  Partial  Fulfillment  of  the  Requirements 
for  the  Degree  of  Master  of  Engineering 

SCATTERING  OF  LIGHT  FROM  NONSPHERICAL 
DIELECTRIC  PARTICLES: 

SPHEROIDAL  CLADOSPORIUMS  (2ym) 

CUBICAL  NaCl  CRYSTALS  (4vim) 

By 

Ismail  Kirmaci 
March  1978 

Chairman:  Markus  Zahn 

Major  Department:  Electrical  Engineering 

A digital  polar  nephalometer  has  been  constructed  to 
measure  the  intensity  and  polarization  of  light  scattered 
from  aerosol  particles  of  controlled  shape.  The  scattering 
at  94  angles  (0°  to  180°  in  1.8°  intervals)  is  digitized  by 
a fast  sainple-and-hold  circuit  and  recorded  on  punched 
paper  tape.  The  light  source  is  a dye  laser  pumped  by  a 
120  kilowatts  nitrogen  laser,  providing  approximately  10 
kilowatt  2 nanosecond  pulses  at  a chosen  wavelength  across 
the  visible  spectrum.  The  device  was  calibrated  with  poly- 
styrene latex  spheres.  Sample  data  for  spheroidal  particles 
(2|im  cladosporiums)  , and  cubical  particles  (4)im  NaCl  par- 
ticles) formed  by  a Berglund-Liu  aerosol  generator,  are  shown. 
It  is  found  that  the  backscatter  from  nonspherical  particles 
differs  significantly  from  that  for  spheres  so  that  a tumbled 
irregular  particle  is  not  equivalent  to  a sphere. 


Chai rman 


CHAPTER  I 
INTRODUCTION 


The  experiment  described  in  this  thesis  was  performed 
to  provide  data  to  test  and  compare  existing  theories  of 
light  scattering  from  irregular  dielectric  particles. 

Several  theories  giving  both  exact  and  approximate  solutions 
have  been  put  forth  independently  in  this  ai*ea  in  recent 
years . 

Ward  [1]  has  an  approximate  solution  for  convex  par- 
ticles using  the  Rayleigh  hypothesis.  Voh  and  Barber  f 2 ] 
have  a nearly  exact  solution  using  Waterman's  extended  boun- 
dary value  approach  [3] , properly  treating  the  region  outside 
the  surface  but  using  an  approximate  extension  between  the 
irregular  surface  and  the  inscribed  sphere.  Eiden  [4]  gave 
an  exact  solution  which  treats  irregular  particles  as  spheres 
with  a wave  number  k which  is  an  arbitrary  function  of  posi- 
tion, which  he  applied  to  objects  with  one  rotational  axis 
of  symmetry.  Herschel  Wiel  extended  his  integral-differen- 
tial representation  of  induced  surface  currents  and  their  re- 
radiated  field  from  conducting  objects  to  dielectric  objects, 
namely  thin  disks  [5].  Ward  gave  a solution  using  a general- 
ized boui.dary  value  approach  [6],  but  without  numerical  re- 
sults. Acquista  (7]  has  developed  Shifrin's  integrodif feren- 
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tial  equation  [8]  into  a potentially  useful  approach. 

These  theories  need  to  be  further  expanded,  applied,  and 
compared  to  determine  the  range  of  validity  of  each.  Our 
role  is  to  provide  further  data  of  high  accuracy  and  high 
resolution  against  which  these  theories  may  be  checked  by 
other  workers.  We  have  calibrated  our  apparatus  with 
spheres,  for  which  the  theory  is  unquestioned,  and  have 
gone  on  to  take  data  for  some  very  simple  shapes.  Close 
attention  has  been  given  to  measurements  near  backscatter 
in  the  range  160-178°  where  little  published  data  exists 
(Pinnick  et  al . [9]  have  data  to  170°),  and  where  the 

scattering  is  strongly  affected  by  shape. 

The  experimental  work  is  done  to  obtain  a relation 
between  the  shape,  size  and  refractive  index  of  the  scat- 
terer  and  the  angular  distribution  of  the  scattered  light. 
Accuracy  is  of  great  importance  to  elucidate  the  effects 
of  the  variables  involved.  The  instrumentation  must  be 
carefully  calibrated.  A polar  nephalometer  is  used  to 
measure  the  angular  distribution  of  scattered  light  in  two 
polarizations.  It  consists  of  a light  source  to  illuminate 
the  particles,  a means  of  collecting  information  by  scanning 
the  scattered  light,  detectors  and  a data  recording  system. 

Many  varieties  of  polar  nephalometers  have  been  con- 
structed and  used  and  are  being  used  for  aerosol  measurements. 
Ours  is  unique  in  its  angular  range,  digital  data  acquisition 
system,  and  tuneable  laser  source.  The  information  obtained 
is  in  form  of  intensity  and  polarization  of  the  scattered 


light  which  later  is  used  to  find  the  properties  of  the 
scattorer.  Accuracy  and  economy  were  the  principal  con- 
siderations in  designing  the  experimental  apparatus. 


The  history  of  polar  nephalometers  goes  back  to  Elter- 
man's  work  (1952)  with  searchlights  which  could  be  termed 
an  early  "polar  nephalometer. " In  1960  B.  S.  Pritchard 
and  W.  G.  Elliott  [10]  used  tungsten  filament  lamp  with  a 
calcite  polarizer  as  a source  with  a synchronous  motor  to 
take  the  scan.  Angular  resolution  was  said  to  be  2°  to  7° 
chosen  according  to  sensitivity  measurements.  Angle  range 
was  not  given.  One  photomultiplier  was  used  for  detection. 
Information  was  recorded  on  a chart  recorder.  Calibration 
was  done  with  a then  new  procedure  which  employed  a diffusing 
plastic  screen  of  known  reflectance  and  admittance  as  a 
standard.  In  this  case  the  sample  volume  was  outside  in 
the  atmosphere. 

The  formulation  modeling  of  light  scattering  pro- 

perties obtained  from  experimental  data  was  tried  by  Kurt 
Bullrich  [11]  in  1964.  He  was  interested  in  the  scattering 
function  and  polarization  of  the  light  scattered  by  a volume 
of  air.  His  light  sources  were  a 500  watt  high-pressure 
mercury  lamp  for  short  waves  and  a 1000  watt  high-pressure 
xenon  lamp  for  visible  and  near  infrared  regions.  He  used 
a "Waldram  iris"  to  keep  the  sample  volume  constant  and 
eliminate  the  background  illumination.  He  also  used  a photo- 
multiplier for  detection. 
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Reiner  Eiden  [12]  studied  the  polarization  (19G6)  to 
determine  the  imaginary  part  of  the  refractive  index  and 
aerosol  size  distribution.  He  used  a polarized  xenon  pres- 
sure lamp  with  uniform  spectral  emittance  as  a source  to 
scan  an  angle  range  from  50°  to  160°  in  10° intervals . A 
photomultiplier  tube  was  used  for  detection.  He  measured 
the  aerosol  size  distribution  for  the  particle  radii  from 
. 1 to  3 ym,  and  pointed  out  that  70%  relative  humidity 

was  a dividing  line  for  aerosol  growth  because  of  the  con- 
densation of  water  vapor  upon  the  surface  of  the  aerosol  i 

particles.  Using  three  different  wavelengths  (.443,  .548, 

.639  micrometers)  he  stated  that,  the  accuracy  of  the  mag- 
nitude of  the  refractive  index  was  1%.  His  results  were 
affected  by  the  turbulence  of  dust  in  the  air  since  he  worked 
in  the  actual  atmospheric  conditions. 

The  first  careful  laboratory  experiments  on  irregular 
particles  were  carried  out  by  A.  C.  Holland  and  Gagne  [13] . 

They  used  a tungsten  halogen  lamp  and  a recording  polar 
nephalometer  consisting  of  two  identical  telescopes,  one 
fixed  at  11°,  the  other  scanning  a range  from  18°  to  166°. 

The  particles  were  carried  by  a jet  of  dry  nitrogen  gas. 

The  calibration  procedure  was  based  on  using  a standard 
sheet  of  known  reflectivity.  Photomultipliers  were  used 
for  detection.  The  particles  used  were  5ym  silica  with  a 
refractive  index  of  1.55.  Measurements  were  made  at  two 
wavelengths.  They  concluded  that  the  light  scattering  for 
a polydispersion  of  randomly  oriented  irregular  particles 
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was  very  similar  to  that:  for  spherical  particles  except 
near  back-scat tor ing  where  the  differences  wore  striking. 

Studies  collecting  more  information  on  complex  index 
of  refraction  of  atmospheric  aerosols  have  been  made  by 
G.  Grams  [14]  from  1973  to  date.  His  polar  nephalometer 
moves  the  chopped  He-Ne  laser  light  source.  The  angular 
range  is  15°-1(>5°.  Angular  resolution  is  limited  to  .5°, 
by  photomultiplier  field  of  view.  In  these  airborne 
studies,  particles  are  also  collected  with  am  impaction 
device  for  microscopic  analysis.  Mie  theory  is  used  to 
fit  curves  for  various  complex  indices  of  refraction  and 
size  distribution  functions. 

H.  H.  Blau  et  al . |Jr>]  constructed  a polar  nephalo- 

meter in  1970  to  study  spherical  Mie  scattering.  A He-Ne 
laser  was  used  to  scan  an  angular  range  of  10o-170°.  A 
photomultiplier  was  used  for  detection.  The  seal terers 
were  suspended  in  an  electric  field;  size  ranged  from  7.5pm 
to  11pm  in  diameter  and  indices  of  refraction  from  1.55  to 
l.f>3.  Calibration  procedures  were  not  mentioned  in  this 
paper.  Blau  et  al . [16]  constructed  another  nephalometer 

to  determine  water  droplet  concentrations  and  size  distribu- 
tions. In  this  case  a Ga-As  diode  laser  with  . 6W  Cw  output 
at  8500A0  was  used,  the  scattering  was  measured  at  30°,  45° 
and  60°.  Calibration  was  made  by  a known  size  droplet  from 
a hypodermic  injection  needle.  A photomultiplier  tube  was 
used;  the  pulse  height  was  taken  as  proportional  to  projected 
area,  and  pulse  rate  as  proportional  to  particle  density. 

In  both  cases  the  angular  resolution  was  less  than  1°. 
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A very  fast  instrument  for  measurement  of  light  scat- 
tering was  made  by  Gucker  et  al . [17] . The  system  used  a 

O 

He-Ne  gas  laser  with  50MW  output  power  at  6328A.  Scanning 
was  made  from  7°  to  173°  and  187°  to  353°  in  14.5  msec. 
Photomultipliers  were  used  for  detection.  The  instrument 
was  calibrated  with  polystyrene  latex  particles.  More 
complex  instruments  have  also  been  made,  for  example  that 
of  Arlon  J.  Hunt  and  Donald  R.  Huffman  [18] . They  used  a 
piezo-birefringence  modulator  to  modulate  the  polarization 
state  of  the  incident  light  beam,  which  was  a collimated 
high  pressure  Hg  lamp.  Detection  was  done  by  a photomul- 
tiplier tube  covering  an  angle  range  of  0°  to  168°.  Their 
results  for  latex  spheres  agreed  with  Mie  theory.  They 
concluded  that  the  fluctuations  were  due  to  the  additional 
particles  in  the  latex  sols  that  were  not  specified  by  the 
manufacturer. 

Many  different  forms  of  aerosol  sizing  and  counting 
instruments  are  commercially  available.  The  problems  with 
these  instruments  are  calibration  and  maintenance.  With  the 
utmost  care  still  the  accuracy  of  10%  in  sizing  becomes  an 
optimistic  estimate.  While  this  continues  to  be  a practical 
problem  that  has  not  been  solved  yet,  there  is  a reasonable 
agreement  between  the  integrated  scattering  coefficient  for 
irregular  particles  and  equivalent  spheres  fitting  the  same 
size  distribution.  So-called  integrating  nophalometers 
measure  all  the  light  scattered  ray  between  60°  and  120°. 
Irregularity  of  the  scattering  particle  has  little  effect  on 
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the  total  scattering  into  this  angular  range.  A comparison 
of  this  is  given  in  a recent  paper  by  R.  G.  Pinnick,  D.  E. 
Carroll  and  [).  J.  Hofmann  119] . The  objectives  of  this 
experiment  are  similar,  but  the  experimental  techniques  are 
different.  Pinnick  et  al.  used  an  argon  ion  laser  of  2W 
output  power  at  514. 5nm.  All  measurements  were  made  at 
this  wavelength.  The  scattering  was  collected  throughout 
an  angle  range  of  10°  to  170°.  Several  monodispersions  of 
polystyrene  latex,  salt,  and  sodium  chloride  were  sampled. 
They  concluded  that  the  scattering  might  be  modeled  with 
some  success  using  Mio  theory  for  spherical  solid  particles 
of  equal  cross-sectional  area,  providing  the  particle  size 
parameter  is  less  than  about  5.  For  larger  particles,  they 
found  that  Mie  theory  overestimated  the  scattering  in  non- 
forward angles. 

This  particular  research  aims  to  help  find  a better 
model  for  light  scattering  from  irregular  shapes  by  exper- 
imenting with  r- variety  of  shapes  and  sizes  of  aerosols  and 
wavelengths  of  incident  light.  At:  this  point  of  progress  we 
were  only  able  to  exoeriment  with  three  shapes  of  one  size 
each  using  one  wavelencrth  of  incident  light. 
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I = 

Q = 
u = 

V = 


E y j and  E| , U 

k is  the  wavenumber  of  incident  light.  The  scattering 
matrix  ( ^ ) is  a 4x4  matrix  which  can  be  reduced  to  less 
than  16  elements  depending  on  symmetry  considerations  in 
various  scattering  situations.  These  elements  can  then  be 
related  directly  to  the  intensity  functions  obtained  from 
Mie  theory.  The  intensity  functions  are  given  in  terms  of 
the  particle  comp lex  refractive  index  m,  size  parameter  x, 
which  is  the  ratio  of  particle  circumference  to  the  wave- 
length and  the;  angle  from  forward  scattering  direction,  0. 

For  -cat*  ring  in  a medium  that  is  isotropic  and  free 
of  optical  activity  the  scattering  matrix  is  found  to  be 
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//  - parallel  component 
[ - perpendicular  component 
* - complex  conjugate 


and  V characterize  ellipticiil  polarization, 
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For  single  homogeneous  spherical  particles, 

S11  = 

S22  = 
s12  " 0 

S33  ~ S 4 4 = li(m,x,0) 

S.,4  — D (m, x,  0 ) 

The  intensity  functions: 

1^  characterizes  the  perpendicular,  I2  characterizes 
the  parallel  polarizations  where  L and  D are  for  elliptical 
polarization. 

This  theory  for  the  characterization  of  polarization 
is  not  shape  dependent.  The  values  of  the  elements  depend 
on  shape.  Measurements  of  the  average  intensity  of  light 
scattered  from  uniform  polystyrene  latex  particles  were  made 
for  comparison  with  Mie  scattering  theory  {see  Appendix  II). 
This  is  used  for  calibration. 

The  scattering  of  electromagnetic  waves  by  an  ellipsoid 
of  arbitrary  size  and  optical  properties  may  be  solved  exactly, 
using  the  method  of  separation  of  variables.  The  procedure 
would  be  to  formulate  the  problem  in  ellipsoidal  coordinates 
and  express  the  solution  of  the  wave  equation  in  a series  of 
ellipsoidal  harmonics.  However,  Stevenson  used  a dif- 

ferent approach  by  showing  that  the  general  solution  of 
electromagnetic  scattering  problems  can  be  expressed  as  a 
power  series  in  the  ratio  of  the  dimension  of  the  scatterer 


to  the  wavelength  (see  Appendix  III).  We  used  this  theory 
to  analyze  scattering  by  cladosporium  spores. 

While  the  exact  solutions  give  a means  of  standard 
to  check  the  experimental  work,  solutions  for  irregular 
particles  are  not  yet  well  known. 

Further  studies  have  been  made  on  the  subject  which 
will  not  be  explained  in  detail  since  it  is  beyond  the 
scope  of  this  thesis. 


CHAPTER  III 


EXPERIMENTAL  APPARATUS 


A polar  nephalometer  is  constructed  to  digitize  and 
record  light  scattering  measurements  (Fig.  1) . It  consists 
of  a scattering  chamber,  nitrogen  pumped  dye  laser,  digital 
data  recording  system  and  aerosol  generator.  The  tumbled 
monedispersions  generated  by  the  vibrating  orifice  aerosol 
generator  are  illuminated  by  the  laser  beam  in  the  scattering 
chamber.  The  scattered  light  is  collected  by  the  photomul- 
tipliers at  several  angles.  The  detected  signals  are  sam- 
pled and  recorded  through  a digital  data  acquisition  system. 

Scattering  Chamber  (Fig.  2)_ 

The  scattering  chamber  is  a light  tight  box  of  dimensions 
( 34x43x78  cm)  made  of  .5  cm  thick  aluminum.  It  is  lined 
with  black  velvet  to  prevent  any  reflections.  The  aerosols 
flow  into  the  box  through  a 7 cm  diameter  hole  and  are  re- 
moved through  two  outlets  which  have  negative  pressure  to 
the  chamber  loccited  10  cm  above  and  10  cm  beneath  the  in- 
let: hole.  The  aerosols  are  illuminated  by  the  incoming  0.5 
cm  dia.  f/12  laser  beam  about  8 cm  from  the  injection  port. 

A 1.5  cm  portion  of  the  7 cm  x 0.5  cm  roughly  cylindrical 
illuminated  volume  is  viewed  by  the  detector  optics  and  is 
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called  the  interaction  region.  The  laser  beam  enters  the 
box  through  a rectangular  slit  of  .5x2.2  cm  focussed  (f/12) 
into  the  interaction  region  and  then  it  is  sent  into  a light 
trap  (Fig.  3) . The  light  scattered  by  the  particles  is 
collected  by  a semicircular  mirror  arrangement  (Fig.  4) 
which  redirects  the  scattered  light  to  an  axial  mirror 
mounted  on  the  shaft  of  a stepping  motor.  Light  passes  out 
though  the  hollow  shaft,  and  is  analyzed  in  a stationary 
Gian  air  prism.  The  light  passing  through  the  prism  is 
decomposed  into  parallel  and  perpendicular  components, 
which  the  intensities  are  separately  detected  by  photomul- 
tipliers. The  secondary  mirror  has  an  angular  resolution 
of  about  1.3°. 

The  volume  of  the  interaction  region  is  individually 
calculated  for  each  laser  beam.  A spherical  lens  with  12 
in.  focal  length  is  used  to  focus  the  nitrogen  lasers'  UV 
output  at  the  interaction  region.  Different  lens  arrange- 
ments are  used  for  the  dye  lasers  either  to  focus  the  beam 
or  to  prevent  divergence  because  of  the  dimensions  of  the 
light  trap.  A light  stopper  is  used  to  prevent  the  direct 
scattering  of  light  from  the  entrance  slit  at  the  backscatter. 
The  mirrors  used  are  front  surface  mirrors,  silver  coated 
by  vacuum  deposition  on  glass  slides  of  1 mm.  thick.  The 
stepping  motor  is  the  Sigma  Instruments,  Inc.  29A-21  Bi- 
polar R/l  Drive  motor  which  has  a torque  of  115  oz-in.  at 
2.4  amps.  Two  microswitches  are  provided  at  the  0° and  1B0° 
of  the  motor  revolution  to  stop  the  operation  in  case  of 


overshoot . 
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Digital  Data 


(Fig.  5) 


The  data  recording  system  is  a very  fast  digitized 
electronics  which  is  capable  of  synchronizing  the  laser 
pulse,  the  stepping  motor  and  the  paper  tape  punch. 

Laser  triggering  starts  with  a manual  push  button  while 
the  mirror  connected  to  the  shaft  of  the  stepping  motor 
is  at  the  forward  scattering  angle,  approximately  2°  from 
the  optical  axis  of  the  laser  beam.  The  scattering  light 
is  collected  by  the  photomultipliers  and  is  multiplied  to 
a voltage  level  detectable  by  a sensitive  instrument. 

This  voltage  is  fed  to  a sample-and-hold  circuit  with  a 
delay  time  of  10ns.  The  threshold  level  for  this  circuit 
is  2.5  mV.  The  samples  voltages  are  digitized  by  means 
of  a 12-bit.  analog  to  digital  converter  in  40iis,  and  re- 
corded on  a punched  paper  tape  with  a BCD  coded  sequence 
count.  The  12-bit  coding  is  shown  in  Fig.  6.  After  the 
information  is  punched  on  the  paper  tape  a command  signal 
triggers  the  laser,  advances  the  motor  and  proceeds  to 
sample-and-hold.  The  main  time-limiting  factor  of  the 
process  is  the  mechanical  paper  tape  punch.  It  takes  about 
.24  seconds  for  the  paper  tape  to  punch  an  information  set 
consisting  of  the  sequence  count,  and  voltages  for  two 
channels.  The  whole  scan  of  94  steps  is  completed  in  23 
seconds.  The  information  on  the  paper  tape  is  decoded  and 
analyzed  by  using  a computer  program  given  in  Appendix  IV. 


The  A/D  converter  and  the  sample  and  hold  amplifier 
used  are  ADC-12QZ  and  SHA1-2A  of  Analog  Devices,  Inc. 
respectively . 

Operation  (Fig.  7 ) 

The  start  box  is  a push  button  control  with  forward 
and  reverse  switches.  Closing  the  switch  in  the  forward 
mode  causes  the  control  unit  to  send  a signal  to  the  relay 
to  release  the  short  circuit  on  the  photomultiplier  inputs 
and  puts  the  S/H  in  the  sample  mode.  At  the  same  time  a 
short  pulse  of  1.7  ws  triggers  the  laser.  The  information 
on  the  photomultipliers  A and  B are  integrated,  partially 
using  the  capacitance  of  the  lines.  After  a very  short 
sampling  the  mode  on  the  S/H  is  set  to  hold  by  the  control 
unit.  The  FET  control  then  chooses  the  input  with  the 
fastest  hold  rate  to  he  digitized.  The  hold  in  one  input 
is  10  times  longer  than  the  other  by  manufacturers  specifi- 
cations. Upon  the  convert  command  from  the  control  unit 
this  signal  is  converted  by  the  A/D  converter. 

The  information  from  the  A/D  converter  is  stored  in 
either  A or  B latches.  After  it  is  stored  the  control  unit 
sets  the  FET  control  to  the  other  signal,  which  is  then 
converted  and  stored  in  the  other  latch  the  same  way.  By 
this  time  the  five-state  counter  is  cleared  and  the  first 
motor  position  is  put  through  gate  1 to  the  output  gate 
and  is  punched  on  paper  tape  upon  a signal  by  the  control 
unit.  During  the  punching  there  is  a busy  signal  that  keeps 


the  counter  blocked.  With  the  completion  of  the  first  punch 
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the  busy  signal  is  released  and  punch  again  command  is 
given  which  through  gate  2 punches  the  first  6 lines  of 
the  information  stored  in  the  A latch.  This  is  repeated 
4 times  until  all  four  6 lines  stored  in  A and  B latches 
are  recorded.  As  soon  as  the  last  bit  of  information  is 
recorded  and  the  restart  is  activated,  the  control  unit 
sends  a signal  to  the  motor  to  proceed  one  step.  The 
operation  up  until  punching  takes  10.4  ms  but  including 
the  limiting  punch  time,  one  set  of  data  is  taken  in  240ms. 
The  operation,  then,  is  repeated  94  times  which  is  part 
of  the  logic  circuitry  in  the  control  unit.  After  the  94th 
step  the  logic  stops  the  operation. 

The  reverse  mode  on  the  start  box  is  to  take  the  motor 
back  to  its  original  position.  During  this,  none  of  the 
other  circuitry  is  activated. 

Nitrogen  Laser 


The  construction  of  the  nitrogen  laser  which  is  used 
to  pump  the  dye  is  from  a paper  by  P.  Schenck  and  H.  Metcalf 
[20] . The  design  incorporates  a bandsaw  blade  as  multiple 
electrode  structure  to  ensure  an  even  transverse  discharge 
in  flowing  nitrogen  and  produces  superradiant  emission  at 
337.1  nm.  The  laser  channel  consists  of  two  plates  of  glass 
supported  on  an  aluminum  base  (Fig.  7a) , which  serves  both 
for  support  and  electrical  ground.  The  high  voltage  edge 
of  the  laser  channel  consists  of  a copper  bus  bar,  bandsaw 
blade  electrode,  and  copper  spacer.  The  electrode  geometry 
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has  been  altered  to  keep  the  discharge  off  the  glass  walls. 
Twenty  590  pF  capacitors  are  mounted  in  parallel,  along  the 
laser  channel  (Ficj.  7b).  The  laser  channels  is  continuously 
pumped  and  the  nitrogen  gas  input  is  leaked  in  at  a constant 
rate.  The  storage  capacitor  is  charged  to  10-15  kV  through 
a resistor.  A Thyratron  is  used  for  switching.  When  the 
switch  close's,  the  positive  side  of  the  storage  capacitor 
is  grounded  and  the  other  side  is  kept  at  a high  negative 
voltage.  Its  charge,  therefore,  flows  on  the  bus  bar  to 
the  dumping  capacitors,  which  provide  the  energy  for  the 
discharge  in  the  nitrogen.  This  results  in  direct  electron 
impact  excitation  of  the  triplet  stages  of  molecular  nitro- 
gen's second  positive  band  and  superradiant  lasing  at  337.1  nra. 
An  Amperex  5C22  Thyratron  is  used  for  switching.  The  best 
measured  power  outlet  of  the  N~  laser  of  40  mm  Hg.  and  12  kV 
is  120  kW  (Fig.  8a, b).  The  pulse  width  is  2 ns  FWHM. 

Aerosol  Generator 

The  generation  of  monodisperse  aerosols  is  achieved  by 
a Berglund-Liu  vibrating  orifice  aerosol  generator  [21] 
machined  by  University  of  Florida  Engineering  Machine  Shop. 

The  generator  is  based  on  the  instability  and  break-up  of 
a cylindrical  liquid  jet.  A cylindrical  liquid  jet  that  is 
unstable  for  controlled  mechanical  disturbances  produces 
equalized  droplets,  which  are  dispersed  and  diluted  to  form 
monodisperse  aerosols.  The  minimum  liquid  velocity  to  form 
a liquid  jet  from  a capillary  tube  is  given  by, 
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where  o is  the  surface  tension,  p is  the  density  of  the 
liquid  and  LK  is  the  diameter  of  liquid  cylinder  jet. 

A practical  formula  for  producing  uniform  droplets  is, 
1.5D.J  < A < 7Dj  where  A is  the  wavelength  of  the  disturb- 
ance. The  disturbance  in  the  system  is  created  by  a 
piezo  electric  ceramic  connected  to  a signal  generator. 

A syringe  pump  delivers  the  liquid  at  a constant  rate. 
This  liquid  jet  then  breaks  up  into  uniform  droplets  at 
the  frequency  of  the  ac  voltage.  The  uniform  droplets 
then  disperse  by  means  of  a turbulent  air  jet  for  the 
particles.  A schematic  diagram  of  the  entire  flow  system 
and  the  dispersion  system  is  given  in  Fig.  9 and  Fig.  10 
respectively.  The  diameter  of  the  droplet  can  be  calcu- 
lated by. 
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where  Q is  the  liquid  flow  rate  and  f is  the  frequency  of 
the  disturbance. 

The  diameter  of  an  aerosol  made  by  dissolving  a non- 
volatile solute  in  a volatile  solvent  is  given  by, 


6QC 

TTf  _ 


1/3 


where  C is  the  volumetric  concentration  of  the  solute  in 
this  solvent.  Correcting  this  value  fot'  impurities  gives 
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Fig.  5 Block  diagram  digital  control  system. 
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Fig.  7a  Cross  sectional  diagram  of  laser  channel. 

1:  Base;  2;  Dumping  capacitor;  3:  Bus  bar; 
4:  Band  saw  blade;  5:  Plate  glass; 

6:  Aluminum  bar  to  centralize  discharge. 
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Fig.  7b  Diagram  of  discharge  electronics  in  nitrogen 
laser. 
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Fig.  8a  Nitrogen  laser  output  through  photo  diode 
oscilloscope  at  40  ns/div  and  40  nV/d.iv. 
FWHM  - f>  ns. 


Fig.  8b  Dye  laser  output  using  Rhodamine  B at  627  n in. 
20  ns/div  and  40  nV/div. 
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Fig.  9 Schematic  diagram  of  Aerosol  generator 
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Fig.  10  Schematic  diagram  of  the  dispersion  system. 

1:  Transducer  (piezoelectric);  2;  25}.jn Platinum 
orifice;  3:  Teflon  O ring;  4:  Housing  1 grounded; 
5:  Housing  2 live;  6:  Slip  cover. 


CHAPTER  IV 


EXPERIMENTAL  RESULTS  AND  CONCLUSION 

The  data  that  will  be  presented  in  this  chapter 
con’sists  of  experimental  and  theoretical  results  obtained 
from  spherical  polystyrene  latex  particles,  elliptical 
cladosporium  spores  and  cubical  sodium  chloride  crystals. 

The  experimental  apparatus  is  calibrated  by  the  comparison 
of  theoretical  Mie  calculations  for  spheres  and  experimental 
scattering  measurements  from  1.091|jm  polystyrene  latex 
particles  of  Dow  Chemicals  run  LS  1028  E (standard  devia- 
tion . 0082)  . A discussion  of  this  is  in  Appendices  VI  and 
VIII  and  the  statistical  program  developed  is  given  in 
Appendix  IV.  Evaluation  of  experimental  results  and  com- 
parisons with  the  theories  are  fully  automated  by  a series 
of  computer  programs.  The  computer  program  used  for  theo- 
retical Mie  calculations  is  given  in  Appendix  V. 

Several  corrections  are  made  on  the  experimental  data 
because  of  scattering  geometry  and  statistical  error  calcuia 
tions.  In  the  comparison  of  theory  and  experiment,  a set  of 
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two  different,  log  intensity  versus  scattering  angle,  graphs 
are  presented  for  incident  light  plane  polarized  parallel  and 
perpendicular  to  the  plane  of  scattering. 

Due  to  the  fact  that  the  scattered  light  is  collected 
through  180°  by  the  rotation  of  a mirror  and  then  analysed 
by  a stationary  prism,  the  intensity  of  the  light  received  at 
each  photomultiplier  tube  is  a combination  of  the  parallel  and 
perpendicular  components  of  the  polarized  light,  regardless 
of  how  the  incident  beam  is  polarized.  This  phenomena  can  be 
explained  by  a simple  rotation  matrix,  which  can  be  expressed 
as : 


M(0) 


cosG 

j-sinO 


sinO 

cos0_ 


The  analyser  prism  divides  the  scattered  light  into  its 
parallel  and  perpendicular  polarizations;  therefore  when  the 
light  incident  to  the  prism  is  rotated  by  0 we  get: 
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(1 


y yCosO 

■ cosO  - 


+ rl 

V/ 


sinO) 2 
sinO) 2 


Ij  and  I2  are  the  intensities  observed  by  the  first  and  second 
photomultipliers  respectively;  1 j / and  1^  are  the  parallel  and 
perpendicular  polarizations  of  scattered  light.  The  terms  on 
the  right  hand  side  of  the  equations  are  amplitude  squares  that 
arc  directly  proportional  to  the  intensities.  So,  to  recover 
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the  parallel  and  perpendicular  components  of  scattered  light 
wo  have  to  take  square  roots  of  the  observed  intensities  and 
solve  for  I//  and  1^.  By  doing  so  we  get: 

Ij^  - /i^  sinO  + 1 2 cosO 
Jyy  - /T1  cos0  - /T^  sinO 


Because  of  the  involvement  and  complexity  of  the  computer 
program  that  processes  the  raw  data,  instead  of  rotating  the 
obtained  results  we  applied  the  inverse  rotation  matrix  to  the 
theoretical  curves  for  comparison.  Appendix  V shows  the  comput- 
er program  developed  to  perform  this  rotation  on  theoretical 
Mie  calculations. 

The  results  are  then  analyzed  as  1^  and  I-,,  first  and 
second  charnels  respectively.  Several  computer  plotted  graphs 
of  experimental  results  and  theoretical  curves  are  given  for 
each  run  of  parallel  and  perpendicular  polarized  incident  beam, 
(an  analyser  prism  is  used  in  the  entrance  slit  to  achieve 
linear  polarization  of  incident  light) 

The  light  source  used  throughout  this  experiment  is  N2 
laser  pumped  Rhodamine  B (Eastman  # 4453)  which  lases  in  the 
range  of  685 nm  to  616  nm  due  to  its  tuning  and  superradiant 
character.  7\n  oscilloscope  picture  of  this  output  versus 
wavelenght,  using  a monochrometer  and  a slow  sweep  of  the 
scope  by  modulated  z axis  on  time  base,  is  given  in  Fig.  11. 

It  was  impossible  to  detect  the  wavelength  of  the  laser  light 
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each  time  the  experiment  was  run.  The  variables  involved  in 


the  tuning  of  the  dye  laser  and  the  operation  of  laser  caus- 


\ 


ed  the  wavelengths  to  be  slightly  different  at  each  set  of  runs 


Therefore,  we  took  the  best  theoretical  fit  to  the  experimental 


data  as  the  theory  to  be  compared  with.  A scanning  electron 
microscope  picture  of  the  scatterers  in  each  case  are  also 
included  in  the  data  set  to  be  presented. 
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Pig.  11  A picture  of  the  dye  laser  output  versus  wavelength 
taken  through  an  oscilloscope.  The  tail  ends  at  641. Onm  and 
starts  at  640. Onm.  The  first  peak  is  due  to  the  superradiance 
of  Rhodamine  B. 


Corrections  For  Experimental  Data 


There  are  several  factors  in  the  experimental  data 
gathering  that  effected  the  results  directly  and  indirectly. 

It  is  impossible  to  separate  and  isolate  these  problems.  The 
major  correction  performed  on  the  raw  data  takes  into  account 
the  fact  that,  the  scattering  is  from  a volume  of  several 
particles.  In  addition  to  this,  the  particles  are  moving 
perpendicular  to  the  light  path.  The  velocity  of  the  particles 
are  controlled  to  be  kept  at  relatively  insignificant  levels 
and  the  number  of  particles  in  the  interaction  region  is 
assumed  to  be  constant.  There  are  also  coagulated  particles 
and  tumbling  effects.  Most  of  these  problems  are  easily 
eliminated  by  averaging  a large  number  of  runs  for  each  data 
point.  The  raw  data  had  to  be  corrected  for  the  volumetric 
scattering  from  the  particles  in  the  interaction  region. 

This  is  treated  as  an  angular  gain  problem.  The  effective 
cross  sectional  area  is  found  to  be  (Appendix  VI) 

A = 2rL  cos®  RcosQ  + L/2 

sinG  R + (L/2) cosO 

L - length  of  the  interaction  region 

r - radius  of  the  interaction  region 

R - distance  of  particles  from  primary  mirror 


This  is  divided  by  the  maximum  cross  sectional  area  which  is 
2rL  and  multiplied  by  a factor  N (0)  which  is  inversely  pro- 
portional to  the  square  of  the  distance  and  the  fourth  power 
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of  sin  0/2. [26]  As  a result  we  got  an  arbitrary  weighing 
function  which  can  be  expressed  as: 


G (0)  = 

sinO 


RcosO  + L/2 
R + (I./2 ) cosG 


(R2sin40/2) 


i 
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The  observed  values  can  be  expressed  as  a convolution  of  one 
particle  with  G(0)  in  the  real  time.  To  recover  the  data  for 
one  particle,  an  inverse  filtration  process  is  applied  to  the 
data  points  since  the  intensities  we  need  are  already  multipli- 
ed by  G(0)  in  the  frequency  domain.  Fast  fourier  techniques 
are  used  for  this  purpose.  A study  of  this  can  be  explained 
as  follows:  Let  a(0)  be  the  intensities  observed  by  the 
photomultipliers.  This  is  already  a digital  quantity.  If 
i (0 ) is  the  scattered  intensity  of  one  particle,  we  have 

a(0)  - i (0)  * g (0) 

in  the  real  time  domain 
Then  in  frequency  domain: 

A ( 0 ) = 1(0)  G ( 0 ) 

Since 

CO  OO  - i ()f 

/ i ( 0 ) g (t-0)  dO  » / 1(0)  G ( 0 ) e J clO 

—on  —<v» 

to  recover  i(0): 


A 


1(0)  = A (0 ) /G  ( 0 ) 
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i (0)  = / l(0)e  l0td0 

which  is  the  inverse  Fourier  transform  of  1(0). 
Then, 

i (0)  = / - e-1  (0)td0 

G (0 ) 


This  is  found  to  be  the  best  correction  (by  statistical 
analysis)  in  order  to  fit  the  calibration  data  to  theory  and 
is  used  for  all  the  outputs. 

Other  corrections  that  had  to  be  taken  into  account  are 
the  threshold  of  the  analog  to  digital  converter,  and  photo- 
multiplier noise.  The  threshold  for  channel  ] and  2 of  the 
electronic  data  acquisition  system  is  measured  to  be  73mV 
and  G 8 mV  respectively.  An  average  of  0.1  volt  deviation  was 
also  found;  and  the  relative  calibration  of  the  electronic 
system  is  entered  as  a correction  to  the  raw  data. 

The  photomultiplier  tubes  used  in  this  experiment  are 
RCA  type  7265.  Photomultiplier  noise  is  measured  to  be  in 
the  range  of  120mV  to  240mV  into  10  Megaohms,  which  yields 
a current  range  of  . 017pA  to  . 0 2 A which  is  well  below  the 
specified  dark  current.  Comparing  these  values  with  the 
scattering  expectations  (Appendix  VII)  wo  made  sure  that  the 
lowest  detectable  intensity  is  greater  than  the  background 
intensities. 
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The  calibration  for  the  relative  gain  of  the  photomulti- 
pliers are  done  by  dark  current  measurements  and  by  using 
several  optical  fibers.  Actually  this  is  not  needed  in  our 
calculations  since  we  give  the  results  in  form  of  channel  1 
and  channel  2 which  are  the  photomultipliers  outputs.  The 
operational  voltages  of  the  photomultipliers  chosen  are  based 
on  this  criteria  i»nd  the  linearity  of  their  gain  in  the 
desired  region. 

Fluctuations  in  the  laser  output  power  are  not  included 
in  these  calculations  since  the  negligible  variations  were 
corrected  by  a large  sum  averaging  at  each  data  point.  A 
statistical  analysis  of  the  averaging  and  error  calculations 
are  given  in  the  proceeding  sections. 

Experimental  Data  And  Graphs 

The  experimental  data  is  given  by  computer  plotted  graphs. 
The  variables  used  are;  the  size  parameter  x,  which  is 

x = 2irr/A 

where  r-Radius  of  the  scatterer, 

A~  wave length  of  the  incident  light, 
and  refractive  index  R,  of  the  scatterer.  The  theoretical 
graphs  are  chosen  to  be  the  best  fit  to  the  experimental  data. 

The  explanations  of  the  figures  are  as  follows: 

Fig.  12  a,b-The  scanning  electron  microscope  pictures  of  the 
calibration  particle,  polystyrene  latex,  of  Dow  Chemicals  run 
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LS  1028  E.  Particle  diameter  given  by  the  manufacturer  is 
1.091  jim  and  the  standard  deviation  is  0.0082.  The  mean 
particle  diameter  we  determined  from  the  electron  microscope 
pictures  is  1.04  pm. 

Fig.  13  - Theoretical  Mie  calculations  for  Channel  1, 
relative  log  intensity  versus  scattering  angle  with  a size 
parameter  of  5.01  and  refractive  index,  1.592. 

Figures  14  to  17  are  with  scattering  volume  of  5cm  length 
and  .3cm  radius.  Concentration  of  the  particles  yield  a 
number  of  15  for  number  of  particles  in  the  interaction 
region.  Corrections  are  done  accordingly. 

Fig.  14  - Channel  1 experimental  results  for  polystyrene 
latex  of  1.04  pm  diameter.  The  incident  light  plane  is 
polarized  perpendicular  to  the  plane  of  scattering. 

Fig.  15  - Same  as  Fig.  14  before  subtracting  the  background. 
Fig.  16  - Channel  1 experimental  results  for  polystyrene 
latex  of  1.04  pm  diameter.  Incident  light  plane  is  polar- 
ized parallel  to  the  plane  of  scattering. 

Fig.  17  - Same  as  Fig.  16  before  subtracting  the  background. 
Figures  18  to  21  are  with  scattering  volume  of  2cm  length 
and  ,5cm  radius.  Concentration  of  the  particles  are  in- 
creased to  have  approximately  35  particles  in  the  inter- 
action region. 

Fig.  .18  - Channel  .1  experimental  results  for  latex  particles. 
Perpendicular  polarized  incident  light. 

Fig.  19  - Same  as  Fig.  7 with  the  background. 

Fig.  20  - Channel  1 latex  data.  Incident  light  polarized 
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parallel . 

Fig.  21  - Same  as  Fig.  20  with  tho  background. 

Fig.  22  - Theoretical  Mie  calculations  for  channel  2, 
relative  log  intensity  versus  scattering  angle,  with  a size 
parameter  of  5.01  and  refractive  index  1.592  for  latex 
particles . 

Conditions  for  figures  23  to  26  are  the  same  as  those  of 
figures  14  to  17. 

Fig.  23  - Channel  2 experimental  results  for  latex  of 
1.04  urn  diameter.  The  incident  light  plane  is  polarized 
perpendicular  to  the  plane  of  scattering. 

Fig.  24  - Same  as  Fig.  23  before  the  subtraction  of  back- 
ground. 

Fig.  25  - Channel  2 experimental  results  for  latex.  Inci- 
dent light  plane  is  polarized  parallel  to  the  plane  of 
scatterine . 

Fig.  26  - Same  as  Fig.  25  without  subtracting  the  back- 
ground . 

Conditions  for  figures  27  to  30  are  the  same  as  those  for 
figures  18  to  21. 

Fig.  27  - Channel  2 data  for  latex.  Incident  light  is 
polarized  perpendicular. 

Fig.  28  - Same  as  Fig  27  including  the  background. 

Fig.  29  - Channel  2 latex  data  with  incident  light  parallel 
polarized. 

Fig.  30  - Same  as  Fig.  29  before  subtracting  the  background. 
An  example  of  how  these  data  arc  tabulated  is  given  in 
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Appendix  VIII. 

Fig.  31  a,b  - The  scanning  electron  microscope  pictures  of 
prolate  spheroidal,  cladosporium  spores.  Mean  particle 
dimensions  are  determined  from  electron  microscope  pictures, 
to  be  3 ym  x 1.5  ym  x 1.5  ym.  That  is  for  the  ellipsoid 
a = 1.5  ym,  b = .75  ym,  c = .75  ym. 

The  experimental  results  for  cladosporium  are  best 
compared  with  Mie  theory  using  a size  parameter  of  12.1 
and  a real  refractive  index  of  1.12. 

Fig.  32  - Theoretical  Mie  result  for  Channel  1,  relative 
log  intensity  versus  scattering  angle,  using  an  equivalent 
spherical  particle  diameter  of  2.47  ym.  The  wavelength  of 
the  incident  light  is  640  nm. 

Fig.  33  - Channel  1 experimental  results,  for  cladosporium. 
Incident  light  plane  is  polarized  parallel  to  the  plane  of 
scattering. 

Fig.  34  - Same  as  Fig.  33  with  the  backgrounds. 

Fig.  35  - Channel  1 experimental  results  for  cladosporium. 
Incident  light  plane  is  polarized  perpendicular  to  the 
plane  of  scattering. 

Fig.  36  - Same  as  Fig.  35,  including  background  values. 

Fig.  37  - Theoretical  Mie  calculations  for  Channel  2. 

Fig.  38  - Channel  2 cladosporium  data.  Incident  light 
plane  polarized  parallel  to  the  plane  of  scattering. 

Fig.  39  - Same  as  Fig.  38  plus  background. 

Fig.  40  - Channel  2 cladosporium  data.  Incident  light 
plane  is  polarized  perpendicular  to  the  plane  of  scattering. 
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Fig.  49  - Same  as  Fig.  48  with  backgrounds. 

Fig.  50  - Theoretical  Mie  results  for  Channel  2. 

Fig.  51  - Channel  2 experimental  data  for  NaCl.  Incident 
light  parallel  polarized. 

Fig.  52  - Same  as  Fig.  51  with  backgrounds. 

Fig.  53  - Channel  2 experimental  results  for  NaCl.  Inci- 
dent light  perpendicular  polarized. 

Fig.  54  - Same  as  Fig.  53  with  backgrounds. 
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Theoretical  Mie  results  for  channel  1,  x=5.01,  R=1.592.  Relative  Log 
intensity  versus  scattering  angle. 
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Experimental  polystyrene  latex  scattering;  incident  light  perpendicular 
polarized.  Channel  1. 
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Fig.  15  Data  of  Fig.  14  before  subtracting  the  background. 
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Fig.  16  Experimental  data  of  channel  1 for  polysyrene  latex  of  1.04pm  diameter. 
Incident  light  parallel  polarized. 


Fig.  17  Data  of  Fig.  16  before  subtracting  the  background. 
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Latex  data  with  incident  light  perpencidular  polarized 
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Fig.  19  Data  of  Fig.  18  before  subtracting  background. 
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Fig.  20  Latex  data.  Incident  light  parallel  polarized 
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Fig.  21  Data  of  Fig.  20  before  subtracting  the  background. 
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Fig.  22  Theoretical  Mie  results  for  channel  2,  x=5.01,  R=1.592.  Relative  Log 
• intensity  versus  scattering  angle. 
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Fig.  23  Experimental  latex  scattering  data;  incident  light  perpendicular 
polarized.  Channel  2. 


24  Data  of  Fig.  23  before  subtracting  the  background. 


V 


cot  »*• 


w W i br  tfio  'c 

oo  3ciwdvo*£ 

00  '4tf'»»b0t«? 
CO  •»0£0l9lfc-? 
00  209C61W? 

oo  3Ci£s*«&#* 

to  30  Zb^bM?*? 

oo  3CP4/b<,o-* 
oo  •iitsm^i‘2 

oo  3tb^.6<*C-* 
oo  3tc>Wb£lC-? 
00  3Cfet>9*0bM 

w man- 

co  3C9±tuebm 
CO  300bM>9b* 
oo  3C9*i*ic*2 

00  3C0^t>fe6M 
00  309*6146* 
00  30/06416* 
oo  30«;bsct>fc# 
00  3tSb«ti94* 
oo  3Cb00t?S* 

00  3C9l09b£* 
oo  3C/4UCI.CM 
10-  3t*99Cp£*/ 
10-30  IbP*96*t7 
10-304  * » Of  t-tr 

1o-3c^2SU5*S 

oo  jtqfcbbSCM 
00  3tPcVb *?*! 
CO  3CPbb9'bM 
00  ?0lWb9<iM 
CO  3C  1X1  W«#I - f 
00  3C  / b 1 C Wfe  * l 
00  3CCfa94lC*S 
00  3C?b6£bC*2 
CO  3CObO£9t,-I 
CO  3Cf7V.Qt?S6"  I 

88  &«««•■{ 
00  3CPt?<*l44*J 
00  3006 U 69*1 
00  30bb‘,lpSM 
00  3C94Vbb1>#J 
0 0 3CS'0*b?C*t 
00  3C?SJiJi*l 
00  3l^V?9*l£*t 
CO  Hblivop*} 
CO  3CWlb0b*» 


00  30bb‘,£pSM 
00  3C9lVbbD#  * 

0 0 30Sbc'b«?DM 
00  3C?SVlJ£*l 
00  3l‘Cb9*l£M 
CO  H*UV9»*1 
CO  3CWUob*» 

CO  3tt»SU£^M 
CO  -fO/OV  dS9*  t 
CO  3CQC'»t>fcZ*l 
00  3tf>Z*.9fC*2 
CO  30?l6t£C*2 
00  3Cf>»4bOl*? 

00  3CvObtbl*2 
CO  3 COVVtf,dmZ 
00  iOibbbbVg 

00  3L,*b«,Tl>£  •? 

00  30«»lV4J?-< 

CO  3t i9*crf£*? 

00  366*116?*? 

00  3004091*-? 

00  30C.0V  9?  I *? 

CO  309  l 6Uc*6  * I 
CO  3*640*99*1 
CO  3C904b?fc*t 
lO-30Pl<rV4b’S 
l(»-  3C?c-4Scf  V- 
10-30 U / bO  l * 4 
CO  30’t,CltVC*I 

co  30i.bPbe*M 

10-300  l l bib* b 
1 0 - 3 0 l i V b 0 4 * £ 
*0-30^9  U9£  •?• 
CO  3l9&/*»**» 
co  30?bObSflM 

CO  30646460*? 

CO  30?S9lt:£-2 
00  !0<b4b‘.£*? 

CO  •i0o,/C9T/!}** 

00  *-*0?bt’ bbfc*  I 
J0-3C  9**960*9 
00  30  l * 4 * h 0 * ? 

00  304SS9lf? 

00  30  c*  46t?lt>*£ 

00  30S bM64*£ 

00  JO?Io>«»T*fr 
0 0 i,U;*U 

00  306b»9»9*& 

to  3Ctcm>otj*tr 
oo  3C^cV£«,C*Si 

A11SN31NX 

••••HX-OC*0»M 


I 


Fig.  25  Channel  2 data  for  latex.  Incident  light  parallel  polarized. 
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Fig.  27  Channel  2 latex  data.  Perpendicular  polarized  incident  light. 
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Fig-  28  Data  of  Fig.  27  before  subtracting  the  background. 
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Electron  microscope  picture  of  c.ladosporium  spores 
Magnification  is  8200.  The  spores  were  not  pre- 
pared for  electron  microscopy;  so  they  shrunk 
during  the  coating  process. 
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Channel  1 cladosporum  data.  Incident  light  parallel  polarized. 
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Tie.  34  Data  of  Fig.  33  before  subtracting  the  background 
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Fic.  35  Channel  1 cladosporiun  data.  Incident  light  perpendicular  polarized 
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Data  of  Fig.  35  before  subtracting  the  background. 


Fig.  37  Theoretical  Mie  results  for  channel  2,  x=12.1,  R=1.12. 
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Fig.  38  Channel  2 cladosporium  data.  Incident  light  parallel  polarized. 
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subtracting  the  background. 
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40  Channel  2 cladosporium  data.  Perpendicular  polarized  incident  light. 
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Data  of  Fig.  40  before  subtracting  the  background. 
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Fig.  43  Channel  2 cladosporiurn  data.  Incident  light  is  noC  polarized 
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Fig.  45  Channel  1 theoretical  Mie  results,  x=10 
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"Fig.  47  Data  of  Fig.  46  before  subtraction  of  the  background. 
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The  major  source  of  error  in  the  measurements  is  the 

ji 

photomultiplier  tubes.  The  main  problem  we  had  was  the 
variations  in  the  photomultiplier  characteristics  after  a 

i 

long  period  of  operation.  In  this  case  a short  duration 
of  high  intensity  light  causes  the  photomultiplier  to  over- 
load, and  its  recovery  might  take  several  weeks.  The  Chan- 
nel 2 photomultiplier  developed  a permanent  noise  problem 
in  the  midst  of  our  experiment.  The  results  can  be  seen  in 
Figures  38  through  41  and  43. 

For  instance  in  Fig.  3,  the  data  points  at  21.5°, 

37.7°,  86.3°,  95.3°,  138.6°,  and  172,8°  are  out  of  range. 

Although  the  scattering  pattern  seems  to  be  somewhat  simi- 
lar to  the  theory  (Fig.  37),  the  intensities  at  27.0°, 

36.0°,  43,2°,  48.6',  77.6°  and  those  closer  to  the  back- 
scattering  range  are  hard  to  explain  without  noise.  Some 
of  these  points  could  be  discarded  as  scattering  from 
doublets  or  triplecs  but  there  is  quite  a lot  of  them  that 
don't  fall  into  the  error  bracket.  These,  we  had  to  in- 
clude in  our  data  presentation  because  there  was  no  reason- 
able way  to  get  rid  of  them. 

Other  than  this,  the  items  for  error  considerations 
are  as  follows: 

1.  Probability  of  multiple  scattering  in  interaction 
region : 

Assume  particle  density  3 x 10®/  f t^ 
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2, 

$ particles  main  beam  will  hit  ur  <b 
r - radius  of  scattering  volume 
- lengtli  of  scattering  volume 


n — concentration  of  particles 

# particles  = 3 x 106/  ft3  x 3.5  3 x 10"G  —3  * 3.16  x --  :<  5 

cm 

= 16.7 


j>  particles  in  cone  from  point  collision  to  mirror: 
, laser  ]/ 

mirror 

h = 2/ 3cm 
0 = 1/36 . 6 rad . 

2r  = hO 

The  conical  volume  from  particle  scattering: 

l 2.  1 w ,3  .,2 

v = — -nr  h = — — li  0 

3 3-1 

v = 55.8  x 10 


c c 

# particles  55.8  x 10~  x 10.59  = 591.48  x 10~  particles 

Due  to  the  irregularities,  taking  h = 1 cm 

-4  3 

v = 8 . 3 x 10  cm 

1 q 

# particles  8.3  x 10-‘  x 10.59  = 8.7  x 10~  particles 
Probability  of  Multiple  scattering  = 

= 8.7  x 10~3/16 . 7 = . 5 x 10~3 

There  is  .05%  probability  of  double  or  multiple 
scattering . 

2.  Probability  of  scattering  from  coagulated  particles: 


07 


Approximately  10%  coagulation  of  latex  and  20s  of 
cladosporium  and  NaCl  was  detected  from  electron  micro- 
scope pictures.  This  is  not  taken  into  account  in  the 
calculations.  Actually  a 10%  coagulation  might  cause 
one  or  two  orders  of  magnitude  increase  in  the  intensi- 
ties at  a specific  scattering  angle.  Hut.  this  problem 
is  avoided  by  putting  an  error  bracket  to  the  data 
points.  This  means  the  intensities  that  might  be  caused 
by  coagulation  of  two  or  more  particles  are  rejected. 

The  values  for  rejection  are  determined  by  Mie  calcula- 
tion. The  averaging  is  done  for  5 runs  at  a time  and 
there  are  at  least  50  of  these  averages.  Assuming 
there  are  at  least  15  particles  in  the  interaction 
region  at  one  time,  one  data  point  is  the  average  over 
3750  particles.  By  rejecting  one  point  out  of  5,  we 
make  an  error  of  20%; 

The  overall  probability  of  error  for  one  data  point  is 

then  .2  x .100/50  = .4% 

3.  False  count  rate  due  to  Gamma  Source: 

There  might  be  a false  count  in  the  photomultiplier  due 

to  the  Kr.86  source  used  to  neutralize  the  aerosol  par- 

3 9 

tides.  This  is  given  as  10  photons/cra  sec  at  1 me- 
ter. The  unshielded  portion  of  photomultipliers  are 

2 3 

approximately  1 cm  . Therefore  false  count  rate  is  .10 

2 2 

photons/cm  sec  x 1 cm  = 1/msec 

Gate  width  of  photomultipliers  = 20  nsec 


-9  -3 

.*.  False  count  rate  = 20  x 10  sec/10  sec 

= 20  x 10"6  = 1/50000 
1 false  count  in  50,000  shots 

4.  False  count  rate  due  to  electronic  command  signals: 

The  digital  data  acquisition  electronics  have  several 
discrepancies  in  recording  the  data.  A high  error  rate  has 
been  detected  especially  in  the  first  and  the  last  data 
points.  This  is  due  to  the  start  and  stop  commands  fed 
into  the  system  during  operation.  The  error  rate  is  exper- 
imentally determined  to  be  as  high  as  80%  at  the  first  and 
last  data  points  which  are  9°  and  176.4°  respectively. 
Within  the  rest  of  the  92  data  points  there  is  a 20%  chance 
of  the  occurence  of  a false  count.  This  happens  when  the 
paper  tape  punches  an  undetermined  code  which  might  be  a 
line  less,  a line  more  or  a noniterated  value  as  the  count. 
In  this  case  that  point  is  regarded  as  an  "out  of  range", 
which  is  indicated  in  the  graphs  by  setting  the  scattering 
angle  off  with  respect  to  the  rest  of  them.  This  data 
point  is  then  discarded.  This,  of  course  causes  an  error 
in  the  overall  averaged  value  because  of  less  averaging. 

5.  Error  encountered  in  subtraction  of  background: 

In  obtaining  the  final  values  the  background  at  each 
set  of  runs  is  subtracted  from  the  averaged  raw  data.  The 
accuracy  of  this  procedure  is  questionable  since  sometimes 
the  background  values  are  higher  than  the  experimental  re- 
sults. background  values  are  checked  against  Raleigh 
scattering,  but  even  in  the  cases  where  there  is  poor 
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agreement,  the  final  values  remain  unchanged.  This  is 
possibly  due  to  the  dust  particles  floating  around  in  the 
scattering  chamber  or  those  that  rise  from  the  velvet 
lining  on  the  walls  because  of  the  turbulance.  Neverthe- 
less we  presented  the  data  before  and  after  the  subtraction 
of  background,  to  avoid  this  confusion. 

6.  Accuracy  of  the  values: 

Accuracy  of  the  detected  values  depend  on  the  sensi- 
tivity of  the  experimental  apparatus  and  the  variables  of 
the  experiment.  Sensitivity  of  the  apparatus  depends  on 
the  angular  resolution,  noise  level  and  threshold  of  the 
electronic  system.  Values  below  .2  volts  cannot  be  detect- 
ed due  to  the  A/D  threshold  and  photomultiplier  noise. 

o 

Angular  resolution  is  limited  to  1 by  an  eyepiece  slit  at 
the  receiving  end  of  the  analyser  prism. 

The  volumetric  scattering  element  explained  in 
Appendix  VT  causes  the  experimental  values  to  be  an  average 
of  the  intensities  from  an  interval  of  scattering  angles. 

As  a result  of  this  any  two  proceeding  data  points  are  not 
independent . 

Using  the  mechanics  of  the  system  the  error  involved 
at  each  data  point  because  of  the  geometry  can  be  calculat- 
ed as  the  ratio  of  the  volume  element  from  one  particle  to 
the  mirror  and  the  volume  covered  by  the  interaction  region 
and  the  mirror.  The  latter  is  a function  of  the  scattering 


angle;  but  we  will  approximate  it  by  taking  the  moan  value. 
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The  angular  dependence  of  the  effective  distance  is 
expressed  as:  (Appendix  VI) 


l)  = L 


cosO  RcosO  + L/2 
sin0  R + I./2cos0 


Then , 


L cosO 
1 ~ 2R  sTnO 

To  take  the  average  we  have  to  integrate  this  over  an 


interval  of  2 it 

D = dO 

av  4ttR  0 sinO 

The  result  of  this  integral  gives  us  an  indeterminate 
form;  so  wo  take  2 (L2/4nR)  1/tanO  as  the  correction  to 
length  I.;  therefore  the  tolerance  on  the  effective  length  is 

* i.2  * 

2nR  tanO 

The  correction  gives  us  an  average  of  2*0,  I.  = 5cm. 

This  gives  us  L - 5 + 2.6 

The  area  considered  is  2rl;  and  r - 1/3  cm 


A - -t  (5  * 2.6) 


which  is 


A = 3.33  ± 1.74 

The  volume;  covered  by  the;  interaction  region  is  the  volume; 
created  by  the  effective  area  and  the  mirror  facing  it. 

The  volume  comes  up  to  be;: 

Vi  - 36.6  (2  ± .86) 

Volume  of  the  cone  covered  by  a particle  and  received  by 
the  mirror  is: 


V = 36.6  x .17 
c 

The  ratio  of  these  volumes  is  then: 

V /V.  = .17/(2  ' .86)  .10  * .05 

c'  1 
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Therefore  the  accuracy  of  a data  point  due  to  the  volume 
effect  is 

± {10  ± 5)% 

We  will  take  the  mean  of  this  as  10%.  There  is  an  addition 
al  5%  error  due  to  the  variation  in  the  number  of  particles 
in  the  scattering  volume.  Since  these  are  additive,  on 
the  average  we  have  + 15%  accuracy  at  each  data  point.  As 
we  mentioned  earlier  this  is  also  a function  of  the  scatter 
ing  angle,  and  this  can  be  taken  as  tin  additional  + 5%, 
which  diminishes  at  90°.  This  means  there  is  a 10%  error 
in  the  forwardmost  scattering  angles  and  a 4-  20%  error  in 
the  backward  most  angles. 

Conclusion 

The  formulation  and  modelling  of  light  scattering 
properties  of  nonspherical  particles  is  of  particular  in- 
terest in  work  with  the  atmosphere.  Accurate  experimental 
data  is  of  importance  to  be  able  to  successfully  accom- 
plish this  task. 

The  polar  nephelometer  designed  for  this  particular 
experiment  is  capable  of’  digitizing  the  intensity  of 
scattered  light  at  94  angles  with  1.8°  intervals.  This 
covers  an  angle  range  of  167.4°,  whose  limits  can  be 
chosen  to  start  anywhere  from  0°  to  12.6°  from  the  forward 
scattering  direction.  We  chose  the  limits  to  bo  9.0°  to 

1 % 

176.4  because  of  the  minor  difraction  of  1 ight  from  the 


entrance  slit  and  the  reflection  from  the  light  trap. 

These  could  not  be  avoided  since  the  divergence  of  the 
incident  light  beam  was  considerably  high  and  to  decrease 
this  would  have  caused  a noise  level  above  detectable  sig- 
nal because  of  the  angular  resolution.  Considering  the 
fact  that  the  last  data  point  is  not  reliable  because  of 

the  electronic  command  signals,  we  achieved  taking  accurate 

o 

data  up  to  174.6  . To  the  best  of  our  knowledge  this  is 
the  furthest  scattering  angle  range  achieved  yet.  In  addi- 
tion to  this  the  large  number  of  data  points  taken  in  this 
range  allows  better  comparison  with  theory  on  an  angle  to 
angle  basis. 

The  vibrating  orifice  aerosol  generator  used  in  this 
experiment  generates  monodisperse  aerosols  of  various 
shapes  from  approximately  0.5  pm  to  50  pm  in  diameter,  with 
an  average  geometrical  standard  deviation  of  1.014.  The 
wavelength  of  the  light  source  can  be  varied  from  approx- 
imately 400  nm  to  800  nm  according  to  the  type  and  the  con- 
centration of  the  dye  to  be  pumped.  This  gives  us  a wide 
range  of  particle  size  parameters  to  measure  and  compare 
the  scattering  patterns  with.  Due  to  certain  unavoidable 
limitations,  we  only  had  the  chance  to  use  one  wavelength 
of  incident  light  to  examine  three  known  shapes  of  one  size 
each. 

The  disadvantage  of  the  experimental  apparatus  is,  the 
time  that  it  takes  to  scan  the  entire  angle  range  (23  sec) . 
During  this  Lime  there  is  a certain  change  in  the  number  of 
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particles  and  their  orientations  because  of  the  turbulance 
in  the  scattering  region.  This  has  been  partially  over- 
come by  averaging  a large  sum  of  runs  at  each  data  point. 

The  angular  scattering  patterns  of  polystyrene  latex 
particles  were  found  to  be  in  very  good  agreement  with  the 
Mie  theory  for  the  size  parameter  of  5.01,  although  the 
experimental  results  were  a little  smoother.  (Figs.  13  to 
17  and  22  to  26) . The  smoothness  of  the  curves  can  either 
be  due  to  the  bandwidth  of  the  laser  light  (Fig.  11) , or 
the  irregularities  on  the  surface  of  the  particles  (Figs. 

12  a,b) . The  effects  of  tumbling  of  a larger  number  of 
particles  was  significant  at  forward  scattering  angles  when 
the  incident  light  plane  was  polarized  perpendicular  to  the 
plane  of  scattering  (Figs.  18  and  19).  This  was  detected 
at  Channel  2 which  is  the  perpendicular  channel  relative 
to  the  body  frame  of  reference.  This  convinced  us  to  use 
a lesser  concentration  of  scatterers. 

We  found  that  angular  scattering  properties  of  prolate 
spheroidal  cladosporium  spores  and  cubical  sodium  chloride 
crystals  are  not  very  dependent  on  incident  light  polariza- 
tion for  the  size  parameters  12.1  and  10.2  that  were  examin- 
ed. Since  there  are  no  exact  theories  for  these  particles, 
we  compared  them  with  Mie  theory  for  equivalent  spheres. 

The  scattering  pattern  for  the  spheroid  is  much  smoother 
than  the  Mie  results  (Figs.  32  to  43),  although  the  figures 
are  not  very  clear  because  of  the  noise  that  appeared  in 
the  system  prior  to  these  sets  of  runs.  By  following  the 


general  pattern,  we  concluded  that  the  backscatter  pattern 
for  a spheroid  is  significantly  different  than  that  of  an 
equivalent  sphere. 

In  the  case  of  scattering  from  cubical  sodium  chloride 
particles  for  a size  parameter  of  10.2,  the  scattering 
properties  were  found  to  be  smoother  than  those  predicted 
for  an  equivalent  sphere  by  Mie  theory.  The  patterns  at 
forward  scattering  angles  are  in  good  agreement,  but  the 
nonforward  angles  are  overestimated  by  Mie  theory  (Figs.  45 
to  54) . The  difference  in  backscattering  patterns  with 
respect  to  equivalent  spheres  are  much  greater  for  the  pro- 
late spheroid  than  the  cubical  particles.  The  smoothness 
of  the  experimental  results  can  be  explained  by  the  varia- 
tions in  the  laser  output  and  the  concentration  of  scatter- 
ers  since  the  forward  scattering  patterns  are  in  good  agree- 
ment with  the  theory. 

In  general,  we  conclude  that  the  backscatter  from  non- 
spherical  particles  differs  significantly  from  that  for 
spheres,  so  that  a tumbled  irregular  particle  is  not  equiv- 
alent to  a sphere. 

We  refrain  from  further  comment  on  the  subject  since 
theory  is  not  our  task,  but  it  is  necessary  to  mention  that 
the  results  of  this  experiment  are  important  for  the 
interpretation  of  measurements  of  atmospheric  aerosols  and 
for  successful  modelling  of  scattering  from  tumbled  non- 
spherical  particles. 


ELECT ROMAGNET I C SCATTERING 


The  discussion  of  light  scattering  theory  starts 
with  Maxwell's  equations  [22],  which  form  the  basis  for 
all  electromagnetic  phenomena.  The  set  of  four  equations. 


V‘D  = p 


3D 

VXU  = J + - 


3B 

V • B = 0 VxE  + -x-  = 0 


are  known  as  Maxwell's  equations,  where  D = dielectric  dis- 
placement; B --  magnetic  induction;  H = magnetic  field  in- 
tensity; E = electric  field  intensity;  J = current  density; 

P = charge  density. 

For  macroscopic  nedia  the  dynamical  response  of  the 
aggregates  of  atoms  is  described  by  the  constitutive  rela- 
tions which  connect  D and  J with  E and  H with  B,  e.g., 

D = e!i*  >1  = a£»  K = pH  for  an  isotropic,  permeable  conducting 
dielectric  where  c = dielectric  constant;  o = conductivity; 
p = permeability. 

The  nonhomogeneous  equations  for  each  component  of  the 
damped  wave  motion  then  become  [20], 


w2..  ~ an  32h  n 

V H - op  - ey  — = 0 

9t  3t2 


However  in  nondissipative  media  the  homogeneous 
undamped  wave  equations  are  expressed  as  (a=0)  , 


V2E  - evt  M*-  = 0 
3t 


V2H  - eh  Mr  = 0 
3tz 


The  velocity  of  light,  in  veirious  media  is  given  by 

-1/2 

corresponding  values  of  (eh)  » since  light  itself  is  a 
form  of  electromagnetic  radiation. 

For  a harmonic  monochromatic  light  wave  in  n homo- 
geneous medium  the  vectors  E and  H are  perpendicular  to 
each  other  and  to  the  direction  of  propagation.  The  scalar 
components  are  related  by 


E -E 

-2S  = — * = Z = 
Hy  Hx  ° 


' 

:-io/w 


where  ZQ  is  the  intrinsic  impedance  of  the  medium  for  plane 
waves;  i is  /-T.  The  associated  electrostatic  and  magnetic 
energies  are  given  by 


U + T = | | ( E | E | 2 + H|H|2)dV 


! 

) 


I 
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For  a periodic  field  such  as  the  electromagnetic 
wave  there  is  a flux  of  energy  crossing  a unit  area  per 
unit  time  which  is  given  by  Poynting' s vector. 


S = ExH 


By  taking  the  time  average  of  this  vector  and 
utilizing  E^/H^ , the  intensity  or  the  flow  of  energy 
crossing  a unit  area  along  the  direction  of  propagation  is 


I = (2Zo)_1  | E | 2 


The  function  u satisfying  equations ; nonhomogeneous  scalar 

wave  eauation  in  dissipative  media  (a>0) : 

V2u  - op  — - ey  ~ 0 

3t  3 1 

and  homoaeneous  wave  equation  in  dielectric  media  (a=0) : 

„2  3 2u  _ „ 

V u - ey  — k ~ 0 

3t 

is  one  of  the  components  of  E or  H.  Then, 


Ex  = Ae 


i (wt-kz) 


Taking  only  the  real  part  of  the  exponential  we  have 


E = Acos(mt-kz) 


where  w = 2iry. 

The  propagation  constant  k is  given  by  the  relation 


.2  2 

k - yew 


lyaw 
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i 


which  can  be  written  as 


where 


and 


“i  ■ “M( 

ei  ' "{*r[( 


1 + 


a2  } 1/2 

2 2 
e co  J 


1/2 


+ 1]} 


1/2 


1 + 


2 2 
e co 


->]} 


1/2 


For  an  undamped  wave  in  a nonconducting  medium  (o=0)  k is 
real  and 


k = cc^  = co(ye)^^  = co/v  = 2tt/X 


where  v is  the  velocity  and  X is  the  wavelength  in  the 
medium. 

For  the  general  case  when  k is  complex,  a complex 
refractive  index  m can  be  defined  by 

E = Aexp[i(iot-k  mz)  ] 

X o 


where  k = 2tt/X  and  m = k/k  = m - im.  . 

o o or  i 

When  there  is  an  absorption,  the  absorption  coefficient 
is  defined  by 


I = Ioexp[-a’zI 

where  the  intensity  of  radiation  decreases  from  .1  to  I 

o 

' over  the  path  length  z. 


r 
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are. 


a'  = knr 


The  components  of  E where  there  is  no  attenuation 

= Aexp [i (t+6^) ] = Acos(t+S^) 

Ey  = Bexp[i (T+62) ] = Acos(t+62) 

by  taking  only  the  real  parts. 

The  locus  of  E is  found  by  eliminating  t from  these 
equations , 

o 2 

E E 2 

- 2 -~£  cos5  = sin^6 

which  is  the  equation  of  an  ellipse, 

B./A  = tana 

where  a is  the  semiside  of  the  circumscribing  rectangle, 
K1 

I = tanx 
Ai 

where  x is  the  semiaxis  of  the  ellipse. 

The  state  of  polarization  is  defined  either  by  a and 
5 or  x and  iji.  The  intensity  is  given  by, 


/E  >> 

rE 

-x-|  + 

g 

Ia  j 

B 

1 = w; 


1 2 2 
2^  (Ai+B?) 


When  the  polarization  ellipse  reduces  to  a straight  line 
<S  = j n (j  = 0,*l,  + 2) 


(j  = 0 , ♦ 1 , + 2 ) 
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and 


E 


1Z  = ( - 1 ) m ® 

K ' L)  A 
x 


When  the  ellipse  reduces  to  a circle,  i.e.,  A=B, 


and 


5 = j 1,1,1 


? 2 2 

E + E = A 
x y 


(in  = ±1,±3,±5,  . • .) 


Another  representation  of  the  polarized  light  is  intro- 
duced by  Stokes.  The  Stokes  parameters  are. 


2 2 2 ? 
s = A + B = AT  + B, 
o 1 1 


S.  - A2  - B2  = (A2»B2)cos2if>cos2X 


2 ...2, 


s._,  --  2ABcos5  --  ( Aj +Bj ) sin2iJ)COs2x 


s3  = 2ABsiniS  = (A2+B2)sin2x 


2 _ 2 , 2 , 2 
s = s . + s _ + S . 
o 1 2 3 

When  6 - 0°  or  180°  the  radiation  is  linearly  pola- 
rized, so  s0  = A2+B2,  s^  = A2-B2,  S2  = *2AB  and  r>3  = 0. 

For  circularly  polarized  light  5 = 90°  or  270°  and 

2 

A=B  so,  s = 's^  = 2 A ; s.  = s0  - 0. 
o 3 12 

One  important  property  of  Stokes'  parameters  is  the 
additivity  for  incoherent  polarized  beams, 


101 


s.,  = 


3i 


The  polarized  and 


natural  parts  are  expressed  as 


s 


o 


s 


n 

o 


The  ratio  of  polarized  intensity  to  the  total  is  then 


o 


s 


n 


o 


[(S^)2  + <SP)2  + (sP)2]172 

s 

o 


Using  a matrix  notation  Stokes'  parameters  of  the  out- 
going beam  must  be  a homogeneous  linear  function  of  those 
of  the 


s' 

o 


“n*o  + a12Sl  + ''  1 3 R 2 + » 


a,,s  + as  + as  + a 
21  o 221  232 


si  = a,,s  + a-,0s,  + a-,0s„  + a 
2 31  o 32  1 33  2 


s'  = a.,s  + a.„s.  + a.,s0  + a 
3 41  o 421  432 


which  can  be  applied  to  scattering  as 
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APPENDIX  II 


GENERAL  THEORY  OF  SCATTERING  bY  A SPHERE  (23] 


The  electromagnetic  properties  of  a particle  may  be 
defined  by  the  incident  wave,  which  has  the  field  vectors 
EL  , It  ; the  wave  inside  the  particle,  Er,  Hr  and  the 
scattered  wave,  Es,  IIf,.  These  vectors  satisfy  Maxwell's 
equations  and  the  boundary  conditions  which  are  respectively 
given  as, 

2 9 E i)  E 

V e - on  rr  - f-u  ^ = 0 

d 9t 

and  ~ 


2 3H  3“H 

V H ~ o n rr;  - ey  ? 


boundar y cond i tions : 

<P.2  r Bj ) • n = 0 

where 


and 


B = pH 

( H 2 - Hl)  • n = L 


where  I.  is  the  surface  current  density.  Once  the  solutions 
for  these  are  obtained,  all  electromagnetic  phenomena  related 
to  the  object  are  known. 

The  electric  and  magnet  ic  Hertz  vectors  n , , and  n ^ , are 
def i ned  by  the  relations 


1 0 3 
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II 

rH 

EQ 

neV  x 

STTj/St 

JL" 

II 

VV-TT^ 

- yc31 2TT_1/3t2 

—2  = 

- y e V x 

3iL2/3t 

—2  = 

VV  * 2L2 

2 2 

- ye3jt/3t  - oyu2 

These  vectors  also  satisfy  the  wave  equations, 

2 a2!-!  E. 

v -l  ~ oy  ITT  " eovl  3t2  e~ 

2 

n 3 Tf  ry 

5 12  ‘ Ol,o  ir  ' e|Jo  -^r  “ 

where  P and  M are  the  electric  and  magnetic  polarizations. 
The  scalar  Hertz  potentials  tt j and  it 2 are  derived  as, 


ir 


1 


-Vtt2  = tt2 


The  spherical  components  of  field  vectors  in  terms  of  these 
potentials  are  then. 


3 (rtr.)  2 

+ E2r  - —2-  + * "l  + ° 


E„  = E,  „ + F. 


32(rTr1)  1 3(r7r2) 

+ k. 


G 10  '20  r 3r30  2 rsinO  3 4> 


1 9 (rn^)  1 9 ^rT,2^ 

E<{>  ~ Elc|)  + 1 ’ 2 rsinG  3r3c|>  k2  r 30 
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3 2 ( r it  ~ ) ~ 

“r  = Hlr  + H2r  = 0 + “^2~  + k ri,2 


J a ( r n 1 ) 3 (rir2) 

H0  = H10  + U20  = -kl  rsinO  ”a$““  + r " 3r30~~ 


H . = H,  . + H 


1 3(riTl)  . 1 3 (r2lT2) 

k,  - — — + 


4.  1 4>  24,  ~1  r 30  rsinO  3r D4> 


where  the  propagation  constant  k = -k^  k2 


k^  = iwe  + o 
k2  - iu> 

The  Debye  potentials  are  the  solutions  of  the  homogeneous 
time-dependent  wave  equation, 


V + k2u’  = 0 


V y 


V I iwt 

where  |i  = y e 


In  spherical  coordinates, 


a 2 r II 1 

. 1 

0 

fsin03n) 

+ i s n 

hr2 

r srnO 

20 

1 ao  J 

r2sin2o  a* 

By  separation  of  variables, 


0 


n = k (r)  0 (0)  4»  (<f>) 

+ ['  k2  - aiajlL]  . *(ri  ■ 0 

dr  ‘ L r J 


I 


JL  (sined0(0l]  + [n(nH 
smO  dO  do  J [_ 


G(O)  = 0 

sin^O  -1 


£ + m24>  (cj))  = 0 

d<J> 


where  m = -n, . . .0, . . .+n  n-integer 


The  solutions  for  the  first  equation  are  Bessel  func- 


tions 


*n(kr)  = (nkr/2)  7 Jn+}s(kr) 

Xn(kr)  - -dlkr/2)1/2  Nn+^(kr) 

where  J , (kr)  and  N (kr)  are  half -integral  order  Bessel 
n+*s  t\+>2 

functions.  The  solutions  to  the  second  equation  are  Legendre 
polynomials 

0 = P^m) (cosO) 

For  the  third  equation  solutions  are  sin(nu|>)  and  cos(m^). 

The  general  solution  of  the  scalar  wave  equation  is  then. 


co  n . . 

tt  v V (m) 

rn  = r l l nn 

n=0  m=-n 


“ I l {cn*n(kr)  + dn(.n(«)}{p,(r’<cosf>} 

n=0  m=-n  v ; v j 

|amcos(m(J>)  + L>msin(m<J>)j 
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m = kx/k2  = >mk0/in2ko  = ml/m2 


Considering  the  Debye  potentials  as  the  incident  tt* 


the  scattered  it®,  it®  and  inside  the  particle  ir^,  we  get 


r r 


rirj  = T-  l i"'1  HWir  *„<>'2r)Pn1)(c°s6lc°s+ 
k n=l 


rn 


1 _ 1 y in_1  fo+.-L-  ib  (k_r)P^  (cos9)sin<{> 

2 72  S n (n+1)  n 2 n 


k n=l 


rlTl  = 


_ \ J in_1  — T^-rr-r  a_L(k-r)Py  (cos9)cos<J> 
, 2 L.  n (n+1)  n^n  2 n 

k n=l 


r"2  ' 


- -V  l i"'1  HTH?rr  bnf-n(k2r)Pn1)(cos0lsin'1, 
k n=l 


and 


r7rr  = A-  y in  1 C (k^P*1*  (cos9)cos$ 

r 1 ,2  L,  n (n+1)  nvn  1 n 

k,  n=l 


rTrr  = JL  y in  1 ^n+j-T  d \Jj  (k,  r)  P'x/  (cos9)sin<t> 
r 2 ,2  L,  n (n+1)  nyn  1 n 


(1) 


k~  n=l 


where  Cn(k2r)  are  the  Hankel  functions 

1/2  (2) 

£n(kr)  = ^(kr)  + ixn(kr)  = (Rkr/2)  Hn+,^(kr) 

The  coefficients  an,  bn  are 

ipn(a)ip^(H)  - 3^(3)^  (a) 
an  = tn(a)*l»„(3)  - 


108 


m<|>n  (u)  >|i^  ( 3)  - <|>n  (8)  ij/ (u) 
bn  _ mCn  («)  (8)  ~ ^n(8)t^(a) 


(prime  indicates  derivative  with  respect  to  the  argument.) 
where  a = k£a  = 21Ia/A  = 2ltm2a/Ao 
8 = k^a  = 2IIm^a/Ao  = ma 
Aq  = wavelength  in  vacuum 
A = wavelength  in  the  medium 
m = m^/n^ 

When  the  measured  scattered  field  is  far  from  the  particle, 
i.e.,  k2r  >>  n,  the  solutions  become  easier. 


t'<k2r)  = i 


. (n+1 ) ik2r 


e 


n ~ik2r 
^n  (k2r)  = 1 e 


With  further  simplification 


g _ H0. 
4»  m2 


-ik~r 

Z co  ^ 

-le  . , v 2n+l 

“V  - “»♦  X imsir 


P ^ ^ (cosO)  dP  lcos0) 

n . , n 

a ; — + b — — 

n sin0  n dO 


E0  = 


HO 

ITU 


ie 


-ik2r 


k2r 


COS * l rT7r 


2n+l 


n-1 


n (n+1 ) 


dP(1)  (cosO)  P(1)  (cosO) 

_n_ . n 

an  dO  n sinO 


Denoting 
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v 2n+l  / 

’1  n4j_  n (n+1)  { an1' n 


(cos0 ) + b t (cos0 ) 
n n 


} 


*2  = J=l  HTiSrr  {anTn(cos9)  + bnTin(cos0)| 


where 


(cos0 ) 

*n(cos9>  - iisr~ 


dp  (coso ) 
Tn(cose)  - — — ae 


Then,  by  Poynting's  theorem,  the  intensity  of  scattered 
radiation  polarized  in  the  0 and  i}>  azimuths  is. 


2 2 
I4>  = X’2  2 lSll  sin  > “ 2 2 1lSin  ^ 

v 4 1!  r 4]|  r 


Jo 


4n  r 


2 2 \ **  2 

2 — 2 I S2  I cos  ^ ~ — 2~2  *2cos  ^ 


2 

4ll2r* 


where  i^  and  ij  are  intensity  functions. 

While  the  exact  solutions  give  a means  of  standard 
to  check  the  experimental  work,  solutions  for  irregular 
particles  are  not  yet  well  know. 

Further  studies  have  been  made  on  the  subject  which 
will  not  be  explained  in  detail  since  it  is  beyond  the 
scope  of  this  thesis. 


APPENDIX  III 


SCATTERING  BY  AN  ELLIPSOID  (A  THIRD  APPROXIMATION)  [24] 

Lot  the  equation  of  the  ellipsoid  be 
x2/a2  + y2/b2  + z2/c2  = 1 

where  a,  b,  c are  the  principal  semi-axes.  Suppose  in  the 
first  place  that  a,  b,  c are  all  unequal,  but  this  restric- 
tion can  be  removed  in  the  final  results.  The  dielectric 
constant  and  permeability  of  the  ellipsoid  are  e , p,  respec- 
tively, the  external  medium  being  a vacuum.  Define  the 
incident  wave  by  the  three  sets  of  direction  cosines, 

(&,m,n) : direction  of  propagation, 

(&j,m^,n^):  electric  vector, 

( ^2 »m2 ' n2^ : magnetic  vector. 

Consider  a steady  state  with  time  factor  e 1U)t.  The 
incident  field  is  then  given  by 

_(0)  . ik(£x+my+nz) 

E = U-pm^n^e  1 

„(0)  _ ,0  m . ik  Ux+my+nz) 

where  k = 2tt/X,  X being  the  wavelength  in  vacuum. 

We  expand  the  incident  field  in  powers  of  ik  in  the 

form 

E(°)  = y E(0)  (ik)p 
p=0  P 

H(0)  = \ H(0)  ( ik) P 
p=0  1 
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Ill 


The  first  three  terms  in  the  expansions  are 


E(0) 


E(0) 

-1 


E(0) 

—2 


( ^2 ' ) 

(£irm^,n^)  (Jlx+ray+nz) 

1 2 
(Hx+my+nz) 


with  similar  expressions  for  . _ obtained  by  replacing 

U , L f Z. 

( «. , mx , nx ) by  U2,m2,n2). 

Let  (E,H)  be  the  scattered  field  and  (e/1^  , the 

field  inside  the  ellipsoid.  Then 


E = 7 E (ik)P,  etc. 

p=0  p 


E q , , E 2 1 are  given  by  the  equations 


*0  - vvo 

EU)  - W*1* 

--0  0 

-1  = -1  + VU1 

H-J1*  = G^1^  + vu;[l) 

E2  = f2  + VV2 

isj1*  = fJ1^  + VV^1* 

curl^  = — Eq 

curlG^1^  = -eE^1' 

cur!F2  = 

curlpj1*  = Hit*1* 

divG^  = divGj^ 

= divF2  = divF^^  = 0 

where  G^ , F2  vanish  at  infinity,  V^,  U^,  V2  are  external 
harmonic  functions,  and  , V2^  internal  harmonic 

functions,  which  satisfy  the  boundary  conditions  on  the 
ellipsoid , 


112 


n x V(VQ  - V(Jl)  ) 

n • V(V0  - ^ 

n x V (U^  - u|i} 

n • V(UL  - pU;*l) 

n x V(V2  - V^l) 
n • V(V2  - evj15 


-n  x E^0) 

-n  • ST 

-n  x (C^  - + H-[0)  ) 

-n  * (Gx  - pG^1>  + H.[0)) 
-n  x (F2  - F^)  + E20)) 

-n  * (F2  - eF^l)  + E^0)) 


where  n denotes  a unit  vector  along  the  outward  normal. 
Further 

n • HjdS  = | n • E2dS  = 0 

where  the  integrals  are  taken  over  the  surface  of  the 
ellipsoid . 

From  these  we  obtain  the  results  as 


n • EgdS  = | 


V0  = AlxIa(^  + A2yIb(?)  + A3zIc(^ 

V^1’  = A^’x  + A^’y  + Ajl!z 

<£l>x  = WE)  - A3yIb(E) 

(G<i))x  = - A» >Z) 

Ux  = BXVTV(^)  + B,XV,TV,(5)  + BiyZIbc(C) 

+ B2zxIca(5)  + B3xyIab(5) 

U1(l)  = B(l)Xv  + B'  (l)xv*  + Bj[ 1 ^ yz  + B2l)zx  + B3l)xy 


— 2 = — 2 + ^—2 


! 


- i h1  x2ia(u  + y2ibU)  + z2icU)  - K5> 


“ 2 B1  x 
L- 


2lac(fJ  + y2lbc(f^  + 3z2lcc(^  " 


+ B2xyIab(a  + ClyzIbc(° 


V2  " Vx  + Vy  + AA 


<J>x  = ~ (x2+y2+z2)xIa(C)  - | xl(£)  - ^-(3a2+b2+c2)xIa(FJ 
(Fji}  >x  = -f-  A1°  (y2  + z2)  - jf  (y2-z2)  + c|l)yZ 


V2  = D1xIa(C)  + D2yIfa(S)  + D3zlc(0 


The  y and  z components  are  given  by  cyclic  permutations 
of  x,  y,  z,  a,  b,  c and  of  the  suffixes  1,  2,  3 in  the 
constants  A^ , h^,  A^,  etc. 

The  notation  is  as  follows:  f,  is  one  of  the  ellip- 
soidal coordinates  of  the  point  (x,y,z),  defined  as  the 
greatest  root  of  the  equation 


2 2 2 

0(C)  H -f—  + -f-  + 1 = 0 

az+C  bz+C  c +C 


Further 


= (a2+v)  (b2+v)  (c2+v)0(v) 


v,v'  = roots  of  equation 


-y — + — + ~f — = 0 

a +v  b +v  c +v 


i 
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« ■ c 


R(u)  = Q(a2+u) (b2+u) (c2+u)] 


x.(«  - r — 

a Jt,  (a^+u)  R (u) 


1.(0  - — 5 

aD  (a  +u)  (b^+u)R(u) 


. . (“  du 

Tv(0  = 2 

V (u-v)  R(u) 

I _ f°°  du 
1 " )0 

j = f ...  d-i 

n '0  (a“+u)R(u) 

, , du f 

i 0 (a^+u)  (b'! 4 u)  R (u) 


du 

r (TTT 


T = T (0) 

V V 

A^  = - (e-1)  f (e)  8.^ 

A<L>  = (2/abc)  f 1 (c ) 

!1M  = [(e-1)  Ia  + Jj] 

p _ ini 1 

(vi-l)T  -(l/vabc)  6(v-v  ) 


same  as  n with  v,  v*  interchanged 


g 1 ( vi ) j-  ij-(b2+c2)  (m2n+n2m)  + m2nc2 
+ n2mb2  + [b2Ib  - c2Ic  - £(b2+c2>  Ub-v]Al 
+ ^(b2-c2)Rf>} 


!L  (U)  - [<b- 


l)(b2+o2)Ibc  + 


(y-l)vabc  ' 


,.(i) Ml 


( jj  — 1 ) v ' abc 


|'in2n+n2m  b2_(,2 


gl(p)  \ abc 


2 Xbc^l 


rl.-I  l 2_  2 1 

+ [“SET  + i <Ib+Io>Ibc],'l 

♦ ^,b2^2)Ibo»«>} 


Cj  = (b2-c2)  [ (a2+v>)  B + (a2+v 1 ) B1 ] 


- (2/3)  |i  (b2-c2)  [ (a2+v)B(l)  + (a2+v  ’ ) B’  ( 1 } ] 


= (1/10)  f1  (e)  | (1-e)  «.1  (a2e2+b2m2+c2n2) 

+ r.  (b2mn2-c2nm2)  + [(2-e)I  - ea2Ia  + 4a/bc]A^ 

2 2 2 (i ) 

+ (cl  +21  )B.  + (1/2)  e v1  (d  +c^)  A.' 

d Cl  -L 

.,2  2.  (i) ) 

- cn (b  -c  ) Bx  ' V 


# *■ 

4 * 


Various  relations  connect  the  elliptic  integrals  (and 
T ,),  of  which  we  mention  the  following. 


] 16 


I + I,  + I = 2/abc 
a b c 


*\  + b2lb  + c2lc  " 1 


Iab  - (Ia'V/<b2-a2> 


3Iaa  + Jab  + *ac  = 2/a  bc 


I 

+ b 1 , + cZ 

I =31 

aa 

ab 

ac  a 

, I I, 

_ 1 

a , b 

+ c 1 — ^ 

2 

2,  ,2 , 

2,  I vabc 

'•a  +v  b +v 

c +vJ 

the 

first  three 

terms  in  the  wave  zone  expan 

3 4 

sion,  we  require  the  coefficients  of  1/R  , 1/R  in  Eqr,  F,^r, 
and  the  coefficient  of  1/R3  in  E2R.  This  necessitates  ex- 
panding the  integrals  1(5),  etc.,  in  powers  of  1/R. 


a = x/R,  ft  = y/R,  y = z/R 

so  that  ( a , ft , y ) are  the  direction  cosines  of  the  radius  vec- 
tor to  the  field  point.  Then 


5 = R2  - (a2a2+b2ft2+c2y2)  + 0 (1/R2 ) 

Hence,  expanding  the  integrals  1(5)  » etc.,  in  inverse  powers 
of  5,  we  find 

1(5)  = 2/R  + (1/R3  Jj-  j(a2+b2+c2) 

+ a2a2+b2ft2  + c2y2~I  + 0(1/R5) 
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IaU)  - 2/3R3  + (1/R5)  [>(1/5)  (3a2+b2+C2) 

+ a2ot2  + b2  p2  + c2y2J  + 0(1/R7) 

I_,b(0  = 2/5r5  + °<1/r7) 

, 

1(0  = 2/5R5  + 0 ( 1/R7) 

For  the  remaining  coefficients,  when  expressed  in  terms 
of  surface  harmonics,  we  find, 

s(0,3)  = _ (4/3)  (A1«+A2(l+A3y) 

S^1,4)  = Pja2  + p2P2  + P3y“  + Mj^Y  + q2Ya  + q3aP 

'*■  + r3Y 


where 


S<2'3)  = (|)A1ct|Ji2a2+b2p2+c2Y2-(|)  (3n2+b2+c2)J 

+ (|)A23[a2n2+b232+c2Y2-(J)  (3b2+c2+a2)| 
+ A3y |a2a2+b232+c2Y2- (j)  (3c2+a2+b2)J 

Pj  = (6/5)  ( (b2+v)  (c2+v)B  + (b2+\>' ) (c2+\>' ) 5'  ] 

qx  - (2/5)  ( (b2-c2)AJ  + 3^1 

r.  » (2/75)  ( (9a2-2b2-2c2 ) Aj  + 51^  - 501^1 

P - k>[|  »«"■«]  ♦ **[-  i S«2-3>  4 A 

+ -A 
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For  the  wave  zone  field  we  have 


E 


0 


9P  + 1 3P' 

30  sine  3<J> 


% = "H0  = 


1 3P  SPl 


[sine  3<|>  30 


4 

where,  as  far  as  the  term  in  k , 


ikR 
e 

R 

ikR 
e 

R 


P = k2(Kia  + K26  + K3y)  + k4|_Lia  + L2g  + L3y 
+ M^a2  + + + ^l3Y  + N2Ya  + N3aS 

* 

- (1/30)  (KlCt  + K28  + K3y)  (a2a2+b202+c2Y2)] 

Here  (R,Q,<f>)  are  polar  coordinates  of  the  field  point, 
while  (af(i,y)  are  direction  cosines  of  the  radius  vector 
to  the  field  point,  and  are  therefore  functions  of  o » <J> 
which  depend  on  the  choice  of  the  polar  axis.  The  con- 
stants , Li , . . . are  defined  as  follows, 

K1  = (§)  (e-l)f1(e)£1 

15L1  = f1(e)|(e-l)  £1[(|)  (6a2-b2-c2) 

- (a2£2+b2m2+c2n2) j + e (b2mn2~c2nm2) | 

+ [f1(e))2£1|(e-l)  [ (e-2) I + ea2Ia  - 4a/ (be)] 
+ e2u(b2+c2)/(abc)  J + f1(e)g1(n)|(Ib-Ic) 
[_(\\/2)  (b2+c2)  (mn2+nm2) -mn2b2-nm2c2J 


x 


ppp 


< 
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-cp (b2-c2)  (mn2+nm2) / (abc) | + [ f ^ ( e ) 1 
x 9i  (y)  (Ib“Ic)  Al  C(E-1)  ki  (^) 

+■  eu  (ey-2)  (b2-c2) / (abc)] 

M1  - i55[(E-1>(Iabnnl+Iacn™l‘2Ibl“l) 

2 2 2 1 
-(2a  JU^-b  nur^-c  nr^) 


2c 


abc (£a2b2) 


15N 


1 = ~ (|)  (p-l)f1(Vi)  (b2-C2)£2  + q1(c) 

2 2 2 2 
x l (e/2)  (b  +c  ) (inn^+nm^)  - b mr^  - c nm-j] 

+ f1(u)g1(e)«.2leu(b2-c2)/(abc)  - (y-Dk^e)] 


where 


kl(c)  * b2jb  ” c2lc  " (c/2) (b2+c2) (Ib-Ic) 


Q = 


'WbcVcaVab' 


4e  (e-1) 


abc 


V ^ 

Aa 

J + — j' 
Ja2b2 


4i: 


.2 


2,2  2 , v 2.2, 

abc  ()a  b ) 


■I 


^ ■ ‘Iab-Iac>/(c2-b2> 


0 [R (u) ] 


■ = r 

.*0 


udu  T 2 . 

= I , - c J , etc. 


[R(u)  1 


3 ab 


v 2 2,22 

l a = a + b + c 
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: I ^ 


I,  I , etc.  are  elliptic  integrals. 

cl 

The  scattering  coefficient  s (ratio  of  energy  scattered 
by  ellipsoid  to  energy  incident  on  ellipsoid)  is  found  to  be 


, . 2,2  2 . 2 2 2 , 2 2,2 . -1/2 
s=(xbc  + m c a + n a b ) 

3 

x | (8/3) k4  l ( I K . | 2 + | K • | 2 ) + (16/3) k6 
l j=l  J J 

r ^ 

X Re  7 (K.L*  + R .£*)  - (1/150) 

L j=l  3 3 11 

x |(3a2+b2+c2) ( | Kx | 2 + | Kj | 2 ) 

+ (3b2+c2+a2) ( | K2 | 2 + | K2 | 2 ) 

+ (3c2+a2+b2)  ( j K3  | 2 + |i<3i2))"|J  + 0 (k7) 

(1)  Spheroid 

If  two  of  the  principal  axes  of  the  ellipsoid  are 
equal,  the  elliptic  integrals  become  elementary  integrals. 
We  find  for  the 
prolate  spheroid  (a=b<c) : 


- 


. 2 2, -1/2  c + (c2-a2)1/2 

I = (c  -a  ) log 2~T72 

c - (c  -a  ) i 


while  for  the 

oblate  spheroid  (a=b^c) : 

I = 2(a2-c2)"1/2  cos'1 (c/a) 


I 


P j 


1 ‘ 

4 • 


Exact  solutions  for  spheroids  have  been  found  by  Asano  and 
Y amamo to  [25], 
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(2)  Sphere 

If  a=b=c,  we  find  easily 

I = 2/a,  I = 2/ (3a3) , I = 2/(5a5) 

a dcx 

J = 2/ (la1) , J'  = 4/(35a5) 

There  is  further  no  loss  of  generality  in  taking  the  z 
axis  as  the  direction  of  propagation  and  the  x axis  as  the 
direction  of  the  electric  vector,  so  that 


n-£^=m2=l 


the  other  direction  cosines  vanishing.  We  then  find 


E„  = 


+ kV 


H„ ' {k2a3(M  co“° + {Si 

x |F(e,n)cos9  + F(p,e)  + 6~(|e+3T  cos2q 

+ 6 (2mW  COS0]}  COS'>  ~I 


E,  = 


t *V 


-Ho  - - {k2a3(Sf + {Si 0080 . 

x |F(e,y)  + F(ji,c)cos0  + 6 (2tr  + 3 j"  cos0 

+ 6tSW  oob20]}  Sin'1' 


ikR 

k 

~R 


where 


F(e,p)  = (3/10) (k+2)  ^ (e^-6r+4+c^p) 


I 
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( 3 ) Perfectly  Conducting  Ellipsoid 

Wo  can  deduce  this  case  by  putting  vi=0  and  making 
I e | -►«»  in  our  general  results.  Thus 


K1  = <2/3IaUl 

L1  = (l/15Ia)  (6a2-b2-c2)  = (a2  Ji,2+b2m2+c2n2)  J 

+ b2mn2  - c2nm2|  + (fc^/1512) (I+a2Ia) 

M1  “ 4?^  [Iabnnl  + Iacmml  ' 2Ibc**-l 

+ — - (2a2  U -b2mm  -c2nn  )1 

abc  (£a2b2)  1 1 1 -I 


Iab1bc  + Ibc'Cca  + IcaIab 
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M,  = 


1 Iabnn2+Iacmm2-2lbcU2 


1 45  WbC+WWab 


K2  r*2 

N,  = £,  + 


1/1L 


1 45Ia  1 (b2+c2)lbc-2/(abc) 


x [b2mn2+c2nm2~ (£1/Ia) (b2Ib-c2Ic)] 


The  remaining  constants  K^K^, 
permutation  as  before. 


are  obtained  by  cyclic 





APPENDIX  IV 

COMPUTER  PROGRAM  TO  DECODE  PAPER  TAPE 
AND  AVERAGE  OVER  STATISTICAL  DATA 


V 


n i vcnsicn 
oivcn'-.jcn 

D!  VF  NSILN 
O I CM 

o i v£nsi  on 

D I MENS!  ON 
O t w r. NS  I nr* 
PivPNsm 
n I M N3  I CM 

O I VCNSI CM 
n I *'F.NSI CM 


AVI  COM.  AV?(94) 

Y I (Mil  ,V,?(04  ; , A VC ( 9 A * ' ANSI  ( 94  > 

h:*ai>c«  (?o.no 

m«.ani  i o a > '?ai  ,va;m(  <m  j .vam.^c^ai  .on  c ta  > ,cri 

OUMI  < OA  ) (f>  A ) 

C»T  A f ( t •?••;*  ) • l'iPr>f  ? II  *1  NOIUM 
AIM  ('*4  > . A02C  r*  4 ) 

M nA  ) • Yf-  I ( 94  ) .Y  r ;?  f 94  ) • YD  I <«4  ) , yon  (V*  | , Y A I < 9 I > • Y AP.  ( *34  ) 
C(Mi9l).A(  !(..)*  A A ( 16  1 
I SK  **  { 2 0 ) » VAl  U«r  c?.  *04  I 
A AM  ( 14  ) •AIM'  C *14) 


Ml  At  Mr  am  I , Mr  AMP 
Rf  M K 

CAM  «rR.RFI,K/HO.O/ 
PPF  *0  . 


YM  \ M| | . 

YM AXP-  I I • 

A C = I • Z /*  7 ♦ A 
A Pa?.  ? 


NUf:  = 2 0 


ISx  - I 

M»  r 16 

I G * I 
I f U 1 = 1 

I N D - I 

i no -1 


I - 1 


N = 9 4 

F A C T “ I . n 
FACIl-^l.F/SF*  1 
00  50  1*1  .N 

*C  « 1*1  . 

c.<  I * I >-  I . 

DC)  51  J’?.NF 

G(  .1  • I I — » LC)AT<  I )*M  J-l  ) 

51  CONTINUE 

AMG(I)  =f*ACTAf 'LOATCIIfAP 
AMOK  I I -I  AC  T | *H.OA  Tl  I ) f Af) 

Of  F 1 C I ) -OF F 
orrpcii-.orr 
OfA  I ( I > -0  » 

DUV? < I)  -0  • 

50  CONTINUE 

c DATA  input: 

C FIP5T  CAIA  CAM)  HF  AOS  (M  TH‘  •»  HUNS  IN  A GMOUP  I HAT  Will  tU  AVI  N AGIO 
C SECOND  C A r • D *»«  FC  I F If.  5 TUP.  NU'1I)!.'I  UACKG»tMiK)  31  AMT  II  S«  IC.IIVLl.Y 

C Tf»|f?C  CM'  IS  MCI  « PUM  TO  Ilf.  SKIM  MHO*  FIHST  #•  ( AKOS  SKIMPS  AM.N*f  INi.i.UOFO 

ne  ioc 5 • iooo)  igps 

loon  rnivAT  ( 1 1 • 20 1 cM 

PI.  tc<  5 • | lot  ) INOll 

tool  rn^vATi  Uiii  o 

Rf  AC  C 5*  1 J00)  ISKP 
DO  IJ  1*1* MOP 


* + 
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125 


PfAt  <5,?0>  t Ml  ADER  ( I « J I , J=  1*  00) 

2 o f r.r.vATtnoAi ) 

I J CflN  UNUf 

I I 

I 2 C r M f NUT; 

*’>!!<(  6 »?l  I C KAO»P(  I * J)  . J - 1 tOO) 

?l  \ nr%#Tf m • • flo a | ) 

4 04  5 CCMtNUt 

IF  ( ITK  *CT  «20)  CC)  TO  4050 
IF  ( |GKP(  I5K)  • EQ.  HOT  1 CO  TO  9050 
40*0  CONTINUE 

imx-.  | r.nsi  ino) 

IT  t I^X.FO.O)  CALL  FX  IT 
IF  I IMf.ini  INO)  . f ) GO  TO  *000 
5050  CfMINUH 

CAt  t 5T*  f S<  STAT  .04  • |MX  .ME  AM  , VEAN?  ,VA9!  ,VAR?  , Of  F l ,C»  F2  . 0) 

5052  CONTINUE 

CALL  AN!  . G .W  . A , Y 1 • V 01  • CS Q. RE C SO. RE C I • IER  I ) 

CALL  •iGILStN.Nf  •MFAM2«G«W  .AA.Y2  #Vl>2  .CSO.RI  CSQ.RLCI  , IF.R2  I 
Rt al cs. i y r f ) ois  lcn.pao 

17  71  rnoNA it -jr* o,5  > 

Rc.  At  (9.  1 7721  CAL  1 .CAL  ?.  PART 
1772  FONVAT (?r 1 0.6.E 16.7) 

DO  f. r Js|*N 
P f " ? • I 4 I 5 

A AM  J)~ANO<J>  *P I / l 40 • 

ARC  C J)  = (4 , *PI  5-LF N + C 05 ( AA N( J>  ) > »L £ N* ( COS ( AA N(  J ) ) + * 2 » ) / ( PI +R AC * S IN ( 

1 A 7 N (J)  ) ) 

55  CCNIfNUE 

DO  <.  C J = 1 , M 

WP  I If  (f.  ,109)  ANGf  J 1 • V l ( J1  .V  2i  J 1 ,ML*  AN  I < J ) • MEAN?*  J I . VA  It  I ( J)  , VAM2C  J) 

I . I Sf'TtL.l  . J>  .5  TAT (L.2. J1  .L - I . I MX! 

100  F C * * A T | • • . 1 7 » 7 . ] > 

AIM  I J ) =C('M  ( VAR)  t J)  ) 

AOS  I J 1-COM ( VA»?< J > ) 

Y I J)) 

Y?  ( J)  -C  CMt  V M J)  ) 

AV  I ( J ) CC  M(  Mr  AN  | ( l))M  ART  | J ) f CAL  1 ) 

Av:  ( J)  (TW(ML  c I > 1 ♦ < A Pf  ( J ) *CAl  ? ) 

9 C Ct  Ml  .‘Ail 

I»  < II  «M  ,GT  .0  ) C.U  TO  o| 

C Al  L GR  AI  M | ( ANC.  • Y I . N.  Y VAX  I • RFR.  PI  t . A.  I 
CC  N T | MUL 

II  I II  »»2  ,r.T  .0  ) GO  1 O 92 

C Al  L C.PAPH1  C A KG  .Y?  , N.N  VAX  7 ,r.|  R.llF  I , K ) 

con  r tour 
n mi NUr 

C A l l graph  \ ( A NG  • AM  l ,N,  YMA  X 1 ,Rf  R ,Rf  I , K ) 

CALI  ;J'APl<I  < ANG.  All?,  N,  YMA  XP.RTR  .Rl  l,K) 

CAll  C.RAPMI  ( A N L»  * A V 1 ,N,Y  VAX.  1 ,Wf  W,  '*(  I • K ) 

CALL  C.R  Af'Hl  ( A NG  , A V ? ,N , YMA  X 2 .Pf  R • RF I , K ) 

1 1 < l • IK  .10.1)  CO  TO  5051 

CALL  MATS!  ST  AT, 94. 1 MX. ML  AN!  , Mf{  A N2.VARI  , V AR2  » DU  M 1 .TUMP.  1 ) 

IT  I 1C.CE .2  ) GO  TO  5051 
I G - I G ♦ I 
CO  TO  5C52 
CCMIMUF 
1C  * 1 
I = hl'9 
I CD  • I CO  M 
I T C.  1 r I CJI) 

HIF  O 

cn  in  I 2 

CPN  I |MUF 

C Al  1 ST  at  S ( ST  AT  ,0  4 , I MX  • ML  A N 1 , ME  AN  2 • V AM  I • V AR  2.  DUM  1 • OL’H  2.0) 

DU  «.0  >1  J - 1 .N 
f)FF  ff  J | -Mc  AN  I ( J ) 

C»-‘r.  C J)  -M  ' AN?  (.1  ) 

A 'I  | ( /»  *rOM(  VA  R|  ( J)  ) 

A i*«  < l|*Cr"(VAPi!l  J ) ) 

AVI  (.1  | :CfM(M=  AMI  ( J ) ) 

A VPI.IPC  «’*«<  MF  AN  21  J ) ) 

w r.  I f.  (f,  ,100  1 ANC  ( J I ,Mf  AN1  ( J I • MF  AN  2 0 ) • VAR  1 { J ) . VAM2C  J ) • ( ST  A T < L . 1 • J ) 
• 5 I a T ( | . ? « J ) • L ~ I » l HX) 

CC.N  I INUF 
INO  ! »!0»1 
I *IN  I 

CC  TO  fill 

CCMINUC 

CAt.l  T'APfR  ( VALUE  1 

1 T r 1 I TOT  » 1 
f S * ' I SK  ♦ I 
CO  Tn  4049 
T NO 


9 1 

°2 

lilt 

C 

C 

5051 

5000 

90 9 l 

50?  0 


ruse  i ion  cnMj  \ i 

CC  X A*  | 0 ./•  (40<I6  . I 
PI  HMI 

cvr 

supmuhnc  i>APrmvAuiE » 

0 I -‘  IISI  CN  LOK  2(  64  I 

V I WINS  ION  VALUE  C2.«14  I 

iKiicm*?  lookup  ii'j6).  inuFtiooi 
in  tt  cm  m,R? 

crEKw  /PUF/ItJUF 

da ia  ir  /i  / 

C A T A irST/O/ 

NT  = '14 
KUX ■ 00 

0"l  13  13  I » I ,NU 

V A L U I 1 .11-0. 

VALLf  I 2,1  I = 0. 

1313  CCMlNUt 

It  € I*  ST  .NE  .01  r, n TO  600 
ro  i?  i-i.iso 

I 3 urKimi  1 1- ->o7 
C INIIIAt  I/I  H'UF 

cn  ir  hi.  ioo 

12  IIILf  II  > >0 

no  ioio  11=1. too 
1=11-1 

1 a = x n>  1 1 . i o » 

Ilia  |/ 10 

I N C = I 111  = 16  IMA 

C i no.(  i a«i  61  nn 

10  10  IP  IX  III’  I INI>  *11=1 

C nuitn  A/n  lcokup  taple 
do  tcr.o  II  1 = 1.2 

ii=iii-i 

DC  5050  I I?  = t .2 
I ? = I I ?-  I 
PC  SOSO  113=1.2 
I 1- II 3-  I 
PU  SOSO  114=1.2 
14=114—1 

Of)  SOSO  115=1.2 
15=115-1 
I'C  sr.so  I I 6-1  .2 
16=  I".-  1 

|Npx  It  •32*l2*l6M3t«.l4*4»IS*?tl6*I 
IVU'.  I 6*  32*  I I (.x  I 4*01  | 3*  At  12*241  I 
L r * T I INCH  ) •>  IV  UF 
5050  Ct  N It  NHL 
If  S I = I 

600  fOMINUE 
I KOI  X - O 

c f ind  nn.i  te  s 

C A I.  L *-  I M>  I 1 F I 
I I I NOS  X-  I K'U  X»  I 

c cn  nr  stlp  « 

if  l ir.r  t I if  i .i;r . iso  I cn  to  60’>6 
i not p=l noKiio ( torn  ir  i ♦ i i. 
if < ikoep.e o . inosx ) on  10  20 
c nua  SrC  urn 

c cfeck  t c sir  if  t me  srr«  » is  hefooe  this  point 
c if  i if  ,i  c.i  » or  to  to 
C no  50  12=1.4 

C I T FST  = IF-I 2 

C IDCF-=LOt:aUPll  r.n  ( I TEST  » »l  1 

C I r I INOEI’.Nt  • I NOEX  | CO  TO  50 

C II 'II.COKUIM  IGF.TI  ITF.ST.5IM  I.EQ.INOCX*!)  60  TO  60 

C fic  COM  I Nile 
C GO  f TO 

C 60  CCNIINUE 
C ME  SYNC 

C II  t IT.LT.2I  Oil  TO  T 1 

C VALU  II  .INOEX-I  >=<50O7<» 

c ir  - ir  - ? 

C T l VALl.t  (2.IN0FX-1  I aO'igig 
C ir-lr-2 

C T 2 CO  TO  20 

C Check  n sr;t  if  nil:  sr ip  » is  aftip  this  point 
to  cc  n 1 1 ' iui" 

C PC  f J I 2= I . 10 

OK  5 0 12  = 1.100 

I TE*I= IF  .1 2 

IF  < ICtTI  nrsr  I.CT  .1561  GO  TO  A09« 

I Non*  »l  cokupc  icr  t 1 1 test  I > 1 1 
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If  < ISUk'P.L  T • INI>I*X  *f>0.  INOI  r»  .GT  . lND|.Xf2  ) GO  TO  «5 
IF  ( !(*•  T ( I TEST  *31  • G T *15  > I CC  TO  0094 
II  (I  HOKUM  IGL  M I Tf  STfS)  H ) 1 I GO  TO  VO 

*5  Cl.MlNUF 
PQ  O N II  NU» 
rMM  hi 

«n  rruKAt  (•  i am  in  trouble*) 

*0*.  a ir  s 1 7t  si 

GC  TO  C096 

90  cc  k it  nue 
c RESVNC 

iMigrn*  inpep-i 

If  I INOFPP.LE  • 0)  GC  TC  2020 
on  n\  I / : | N'V  X , I NObPP 
V4(  IF  ( f ,1  /)  rOO  >99 
Vl  VMA*  12*1/1 
l MU  x - IKOr.P 
2C20  CCMIUUC 
I r = 1 1 r s T 

20  CC  M INUE 
C PUll  O II  C NUMIIFIIS 
Of?  I0O  hl«? 

IF  T ir  M 

IF  ( lGr  T(  ir  )*GT.  156)  GO  TO -6006 
I P 1 ( I OCT  I IP  )•!«&)»! 

W 1 =LTK2  ( 1 PI  ) 

If  = IF  U 

IF  < IGr.T  | |F  I *cr  • 156  I GO  TO  6096 
|P2:MOL)I  IGF  T I IF  ) *6A)M 
V2  I OK  2(102) 

100  VAl.Cl  ( I • INDEX  I »R1  *6*  f«2 
IF-  If  *1 

ir  < ik or  x.nc.mo)  go  to  11 
6C<;s  ccMfar 

C A U.  ( | NO  (IT) 

m;  Ttr-M 

6C9e  ccn  1 1 nor 

IF  ( INin  X*  GL  *NO-2  ) G(1  TO  0005 
I a nr  x o 
c ai.i.  i i Nf>(  ) 

CC  TO  II 
r no 

r.UPM'OT  1NE  FlNO(irp) 

If  s|rr 

VC  60  CUM  1 NtIF 

ir  ( IO  I < |)  ) ,Nt •?£5 ) GO  TO  vovo 
I F ( I 61  Mini  GC  TO  vovi 

vovo  crsiiMur 

IF:  |f  ♦ I 

go  to  vono 

90V  I CON  I INUfc 
If  = !•  ♦ 1 
I r r - i f 
Ft  T 0*5 M 
ENC 

r u ac  i ion  i g<  i ) 
iNliCl.PF?  I OOF  (100) 
tor* km  «i  iri(  ' oo) 
too  IV  AC  i -ICC  ( ir.uf  ( 1 1 • UK  i ) ) 

• CCf*»''HSftUFS1  1UF 
DAO  I *>  Vc  / 0/ 

MIA  MO 
AC  A/>  - N U A' • 2 
J - v.-n  ( I *NUM| 

If  ( J.I  O .0  ) J NON 
If  ( J.N‘  *1  ) GO  10  2 0 
IT  4 I *c.O  • I SVfc  > GO  TO  20 
PI  *0(0.10)  ( l ll(  1 I ) * I l -2*  NUM?  ,2  ) 

) 0 F II^a  I ( *0%  1 ) 

w«t  ii^(g.ioio)  muF 
loio  fn.vM(poiM 
70  ic.-t:n/r(jf 

K*  V F s l 

•*r  u»n 
fnc 

SC’I  FOOT  INI.  GM  APMJ  ( K ,v  • N*V  MAX,  «FR,  OF  I . X ) ' 

D 1 Ml  NS  UN  K ( N ) . Y | N ) ,L  1 NF  ( 1 0 1 ) 

PC  At  I INC  . K 

DATA  OOT .STAR. PLANK/*  ••••♦•  *•  •/ 

PO  1 J - 1 • 1 0 I 
1 l 1 A C ( J ) “OPT 


1 


128 

' 


I 

t ! 


? 


I 


»»pnKf».ioi  t rp  #pf  i ,x 

10  A I ( • | • ,10  X,  *M  - • ,1  * .2,  •>|M*.f'4.3,SXi  *K  = • ,F  A.  1 > 

Wf  f ( ft  .1  1 I 

11  run  .-A  I t *0  • ,/X  , • AN(il  r_  • ,7  X,  • IK1FNSI  T Y • » 5X  , * ,01*  •?2X,,0.I»  ,c7X,M  .0*, 

l?n  x • • I u • *vx  . • loo  • > 

HV  I V c f-  ./  1 L IN?. 

6 rcirv/T  it  | 0 , | , l»>  £ 1 ft,  7,  4X.  1 01  A1  > 

2 f ( 10X  .1 01  A1  ) 

DO  « I S I |N 

IF  (APSCYll))  ,GT.  YMAX ) CO  TO  A 

nu  2 j-  i . i oi 

3 l I MM  J ) aML  ANK 

L I M (M  > -DU1 

J a SC.  0*  C Y ( I > / YMAX  ♦ 1 .0JM.fi 
L !*r  | J )?Sf  AP 

*«|IL  t 6 if.  1 K<  1 1 ,Y  t I ) ,L  IME 

CO  TO  5 

4 * p i Tt  I a , 7 1 k ( n • v t i > 

7 FTP  WAT  (r  1 O.S  |E1C),7,5X«  * OUT  OF  PANGEM 

5 Cr^UNUF 

PE  IUM 
Eftn 

FUNCT  K N I GET  ( 1 ) 

r>  i w * k 7.  i r n i »\r  f 1 1 o > 

data  lSV'/O/ 

-JA  1 A It  ST/O/ 

0 A 1 A I MAX/10/ 

If  I ISVF  ,l0  .0  > CO  TO  MO 
IF  ( I • G T « I i I GO  TO  SO 
10*'  =C  t - I GVF  ) f II  ST 
IF(  iro.LI.il  GO  TO  MOMO 

1 C “ I”  I Of  r ( I OP  1 
It  5 1-100 

I S V t - I 

or  »c  <*'>') 

so  ci  nt  riue 

inrtii-r  vr  i 
Ik  : |i  si  ♦ i ;>n 

IF  flK.LT*  max  ' (*0  TO  TO 
ll>  II  l - I ► - I *AX 

I MU  A,  - | MAX  — I DIET 

or  i j - i , i muk 

J.»  r 

I f r r ( j i \ »'n  ( j.i  ) 

51  CFMtK'MF 

on  r 7 j - i , loir  r 

Jj  - J*  I MI.K 
If  =1*  M 

IIH  r ( j j |-  | G(  I F | 
f?  ( r M f MJF 

I L 5 1 I L S T ♦ 1 OP  - I O 1 f F 

Hr  v ini  » c it  st  ) 

I s vr  = i 

GC  TO  C *}<} 

60  CCNTINUF 

ir.r  t - mrr<  |K) 

1 L S f r-  I K 
t S VF-  I 
GC  10 

QO  C r K T f N(J  E 

on  S'»  j - i • I ma  x 
I S T s I-  1 #J 

oo  |r*tf(j)=|G(IST) 

I G r 7 ! P F r ( ) > 

1 r>  V F - I 
! I <M 
If  s | rt 

one  CCATINUE 
les  3 FtjirvATC  e»IS) 

PI  fUf>N 

0040  C f N T I N\J( 

oil  f»  ( ft  • 1 On*))  IGF  T . I , It  ST  ,1  K • I MUK  , ISVI 

» *•  I T f ( ft  « 1 / 0 ) 

1/0  FOOMATC*  wMAT  APF  YOU  PCI  KC.7  I AM  POM  M l NO  MYSF1.F  OUT  •) 

L AH.  Fx  I T 
F NO 

IMG**  E MY  SNIFF!  MT  NC  T PPf  Cf  n p #?  V A TRANSFTR 

SUP TOUT  CNF  ST  a?  sc  S r AT ,N,  f MX, ME AN I , M{  An/, VAP 1 , VAR  2. CM  1 ,CF  F ? . I SK  ) 
r 1 Vf  MSI  ON  »TA1  CS  ,/  ,«»0  ) .me  AM  <Q4  ) . Ml  AN/  ( 94  ) , V AP  l ( O 4 ) , VAP?(  34  ) 

O I Mf.  N S | f?N  fir  I l ( s * I • ( »»  f ■ / I *>  A > 

DlMrN",irN  V Al  UF  C / • f>  4 ) 
r»l  A | Mf  API  | • Iff  AN." 
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\r  i !r.K  , f f)  , I ) CO  TO  50 

f»r  <*  ')  K . |M* 

C AL  L PAPrp  ( VALUf  I 

no  r;  j j-i,m 
s t a i ( k 1 1 • j > = v al  ur  < i . j ) 

r.T  A T < K . ?,J  MV  ALUH  ( J#  J ) 

CON  1 I NUE 
rf  Nil NUF 
50  COM  I NOE 

do  s?  j = i • n 

me  AM  ( J >=0  • 

V{  AN  * c J ) =0  . 

V^I(  j)  =0. 

V t ff  r ( J ) -O  . 

DO  ^ » <a|  .|MX 

ME  AN | ( J I-MIAU  1<  J ) * S T A f ( K t I . J| 

ME  A KC*  ( J ) *:  M f AN  ? I J > »5  T A T ( K • ? . J | 

CL  MI  NUF 

Mr  AS  | ( J J-rMI  AN  1 C J )/rt  r.ATIlMX) 

Ml  AN?  ( J ) - Mt  AN?  ( J I /M  OAT  ( | MX  ) 

DO  « A K - I • I M x 
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VA  ??  ( J)  -(  Ml*  AN?  | J)  ST  AT  IK.?.  J)  ) * •<?  *V  AR?  < J ) 

cc.n i iMur 

VAf.m>  -V  API  ( J)/M  f)AT  ( I MX  ) 
van;  | J)  -VAO?C  J)  /M  IJAT  ( 1 MX  ) 

Mr  am  ( j i Mt  am  i | j ) -urr  i < j) 
mt  an?( j ) -me  an? c J) -err? ( j) 

If  ( I'f  an  1 ( J > .1.  T . C.  ) Mf  ANK  J|-0. 

IT  ( **E  Af4?c  J I.LT.O.)  MCAN?(  J ) * 0 • 

C(  N linue 

IM  III  M 
END 

SUMMHIT  INF  Mr,!  I 5 C N ,M  • Y • • W • A .YF  I T .D  1 1 .CM  ISOM*  Of  DC  SC  . ME  CCMI  • I CP  I 

WGILS  IS  A Gr.MEMAL  LFAST-SGUAnr  riT  ROUT  INF  WHICH  D C T F R M IMPS  THE 
corrr  icifNTs  ai,  a?,  ...  • am  in  tic  tunc  thin 

YF  t 7 -A  I *GI  { Xi  ♦ Af+iiPlX)  * ...  f AM  * CM  (X  ) 

FIT  lu  THE  IN  rvU  r A»r  AY  Y AT  N DATA  POINTS. 

MOLLS  I’H’riWMS  T ME.  SAME  CALC  ML  at  I on  AS  SUHROUT1  Nfl  (ill  F . WITH  IHh 

rnur*iK(.  ado  i r i rus : 

1.  V>r  tt.MT!,  Ar>f  AI.LCwr.D  AND  ARf*.  COMPULSORY.  TOC  ARRAY  iv  IS  USED 

TC  WT  IC.HI  I HI  DATA.  FOR  AM  UNWf  K.H  U D FIT.  W MUST  CONSIST 

FF  All  CNF.  S. 

2 . TMl  H UNO  T I ON  ••  YE  I T ,0!  T .CIH  SOP  ,Rt  UCSO.REIK  HI  APE  CALCULATED  AS 

m scp  i p i)  nti  nw  . 

c.  IS  AN  INPUT  ARRAY  UlvFN  HY  G ( J • K)  -G  J<  X<  K ) ) ; l.t.  G(J.K)  IS 
If- r jll  I'll  he  T I f_)N  c.J  IN  THE  AMOVE  EXPRESSION  FUR  Yf  l T . EVALUATED  AT 
T HT  KTM  grid  POINT  X(K).  (NOT I THAT  X KTSELf  IS  MJT  NEEDED  AS 
AN  fNFUF.I 

rr  Al  Y(  U)  .G(  M.M  . W (N>  .A  <M  } . YF  I T(  N>  .DIF  (N)  tF  D(  6?5  ) 
t NFUT S - - 

n:  m mmf.r  nr  data  points  in  uf  fit. 

m:  M.rnrn  rr  cui  f-r  icifmts  to  of  octcmminco  (m.i  l.n;  m.lg.25) 
y:  f.-n.rcNT  array  ccm  aiming  y -values  or  t nr  data  to  me  fit. 

r,:  (f,M)-tlfMI  Ml  ARRAY  OF  FITTING  FUNCTIONS.  OESCRIHCO  AMOVE.  • 

*:  N-Hl  MINI  Al*.»»  AY  Ilf  INPUT  W 1 F GHT  S • TfHS  ROUTINE  MINIMIZES 

i*«r  function:  sumc  (f  -y  >•  ♦?  > • 
mi  tpl t r -- 

a:  M*f  t r Ml  V T A'  PAY  OF  CMFF  r It  I ENTS 

Yf  I T : M-ri  1'ir-iT  FPPAY  containing  Tut  VALUES  Of  T MF  FUNCTION 

Yf  If*  StMT  A«f.  ) 

Die:  M-ILi  Mf.'JT  ARRAY  * I TM  ERROR  VALUE.  S.  Y - YF  IT. 

FPIGRR:  (Ml  SOUAF-f  --  SUM(  WW(  YF  IT-Y  ) • ♦?> 

Rff’C'  ):  Rf  P*IC!  n col  SflUAI!*-  --  CHI  SOR/ C N-  M ) 

rffoi:  f-riurcn  Cf  i --  sex*  r (pfdcso  f 

IFF:  f RPC  p IN.'If  ATOP-- 

If  P l>:  FK  f MAI  T*  IIMINAT  ION. 

If  r-1  : MtKVll  «rv  MATRIX  PI  TURNED  FROM  SI»1Q.  PPOPAOLE  CAUSE: 

TUfiCTlUNS  C.1  . USC 3 TO  Of  (P  RATT  G APf  NOT 

I |N‘  API  Y I MJfcP  FNDFN  T • 
ir;. r. m too  oic. 

SO  IT  PR  C (RAMS  CALI  Fl»: 

Give:  IHH  SClf  NT  IF  1C  MMIPOUT  INF 
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DIMENSION  STFETpl 10  1 l.SELTRMl A.  loll  .POL  1 < 10  t ) .PPL  21  lOII.UdOll 
D I ME  NS  I ON  PI  Ml  < 101  1 .PI  NT  2(  1 01  > 

O IMF NS  ION  ELI K 101 1 .FLT21 101 I.ELO t 1 101 > .EL  021 10  1) 

PE  ALA  5 X ,RFP ,RF  I.THETDllOOl.QEXT.OSCAT.CTtlROS.F,  LTRMXlA.  10  0.  2) 
PFR=l.«92 

RFI=0.000 

PI *3. 1 « 1 59  2G5 3 5 09793230*6 
PrAOIS.200 > PM.E.D 
200  MlC!'ATUrlO.«) 

on  loo  1=1  ,toi 
ri  ui  ( i i-o. 

EL02I  I 1 =0. 

100  C0N7INLE 

DO  7 N - 1*101 
7 SIMETiX  M =-l  .0*NM  . n 
DO  n N= 1 ,51 
e Tt'ETOlM  =S  THE  TD<N  ) 

00  25  KKsl.11 
RE  APIS.? SO)  P 
250  FORMAT  IF  10. 6) 

V = E API - 1 (P-PM )/E)»*2/2> 

X=PI*D/P 

CALL  013  M t £ <x .RFR.RF I .T HFTD.51 . OEXT ,OSC AT ,CT nnOS.FLTMMX  1 
DO  9 J=I.S1 
DO  9 1 = I . A 

9 SrLTRMl I , J1=ELTRMX( I.J.l  I 
OO  10  J=[.50 
JJ=51-J 
DO  10  1=1.5 

10  SELTRM1 I,S1*J)=ELTPMXU .JJ.2) 

DO  30  *=1,101 
FLT!lK)  = r*SFLTRMM  ,K> 

ELT2(K)  = r*Sr;LTRM(2.KI 
El  Cl  (K)  =EL01  IK)  *El-T  1 <K> 

EL 02 ( K 1 =ELO?<K) FELT  21 K) 

30  CONTINUE 
25  CCNT1NLE 

00  20  K = 1.1  01 

UIK )=S1 PET OIK  I »Pt/l PO . 

POL I IK) = I E LCI  (K)*fOS(U(Kl )*FL02  <K)*SIN(U(K ) > 1 *»2 
Pni.  2 ( K 1 = (EL02(k  )«COS(  U(K))-EL011K)*SIN(U(K))I»*? 

20  CONTINUE 

on  21  K=1,!0I 

P INI!  ( K ) = AL  OC.  10  1 POL  1 (K>  ) 

21  PINT2(*)=ALCG10<PCL2(XI) 

CAl.l.  GRAPH  I 1 STHETP.  PI  NT1  .1  01 .1  O.O.PFR.RFI  .X  1 
CALL  C.R  APH1  ISTHE10.P1NT2.10  1.10.0  .RED  .Rr  l . XI 
S T CP 
FNO 
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5UMH.JUI  INK  CFME  ( X • Rf  R ,RE  I . T ME  T 0 • JX  • Of  XT  • CSC  XT  •CTpROS*  f.L  TH  MX  I 
SlICIWOU  T I NE  row  COMPUTING  T ME  PARAMETERS  Of  THF  1'l.FC  T ROMAGNE  TIC 

computations  i n i>f ’UHiF  pwrcisinN  arithmetic# 

this  sin'inufiwi  computes  tiif.  capital  a function  nv  making  usf.  of 

OnfcN*APO  WFCUPPfN CF  RELATIONSHIP. 

X J SI/!  PARAMETER  C»'  THE  SPMfPt  , ( ? ♦ PI  ♦ RADIOS  OF  THF 

WAVfltKGTH  Cf  IMF  INCIDENT  RAOUTIOM. 

Rf  J REFRACT  I VF  INDEX  OF  THE  MATERIAL  OF  THE  SPMF RF • C OMPLE X 
OUANT I T V%.FC FN t l WFP  - I * PF | ) 

thet;mj>:  ancle  in  orcnri  s oftwein  the  directions  cf  thf  incioent 

and  1 F*t  SCAT  If  REP  RADIATION.  THE  TO  I J ) IS  < OR  * 00.0 

T»  THCTOIJ)  SHLULO  HAPPEN  TO  OF  GREATER  THAN  00.0.  F.NIF.R  KITH 
SUPPl  E.ME.NT  AM  Y VALUE.;  SEE  COMMENTS  BELOW  ON  R IRMX.. 
jx:  TOTAL  NUMOER  OF  THfTO  FOR  WHICH  THE  COMPUTATIONS  AWE 
RfOUlPFO.  JX  SHOUT  0 NOT  I Xf  FF  O 100  UNI  F SS  THE  n|MINS|()NS 
STATEMENTS  AWf  APPROPRI  ATfLY  MODIFIED* 

MAIN  PROGRAM  SMpUtl)  AISO  HAVE  RFAI.+  A 1 ME  TOl  100  I *F  I TRMXl  A.  100.2)  • 

THF  OIF  INI  TICKS  FOR  THE  FOLLOWING  SYMBOLS  CAN  HC  FOUND  IN  • LIGHT 

SCAT  TIRING  f l V SMALL  PAR  T I CL  F $ . M • C • VAN  OF  MUIST.  JOHN  WILLY  f,  SON 

• INC..  N'.  W VCUK,  I SS  T • . 

orxi:  Iff  I FC  I T NC  Y I AC  TOR  FOW  rxTl  NC  T ION  .VAN  OF  HUt.G  T#l**  14  f.  I2f. 

QSCAT:  EEflfCIEKCV  FACTOR  TOR  SCAT  f F R INC..  VAN  OF  HUt.  S I ,P  . I A C I 2 7. 

cumos:  AVERAGE  I C CS  INC  T NET  A>  * OSC.AT.VAN  OF  » H»L  S I , P . 1 ? R . 

■ELIRMX  ( I • J»K  ) : LLFMfNTS  OF  IMF  1 R AN  5*  'OR  MA  T 1 ON  MATRIX  » .VAN  OF  HUL 

ST.P.3A.A5  L ITS.  I t:  ELFMENT  M SUB  ?..t  ?:  El  I ME NT  M SUB 

I * 3:  a»  m^nt  s sun  ?i*.  i * a:  element  o sun  ?!.. 

FI  Tl  mk  < I , J . I | RflMfl  ‘INIS  THF  ITNFlfM'NT  OF  |Mf*  MATRIX  TOR 
TH!i  ANC.Ll'  TEE  TIM  J)*#  F L T WMX  ( I • J .2  I REPRESENT  S T Ml  Mu  III  Ml  NT 
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APPENDIX  VI 


VOLUMETRIC  CORRECTION 


The  figure  given  above  represents  the  scattering 
geometry.  The  direction  of  incident  laser  light  is  + x,  and 
the  direction  of  flow  of  the  aerosol  particles  is  + y. 
Aerosols  enter  the  chamber  through  a circular  opening  of 
diameter  AB.  They  get  hit  by  the  laser  beam  in  the  inter- 
action region  AB.  The  scattered  light  is  then  collected  by 
the  semicircular  arrangement  of  mirrors  A'C'B'.  The  inci- 
dent light  is  focused  at  the  interaction  region.  The  radius 
of  the  light  beam  is  negligible  compared  to  its  length. 
Therefore  the  problem  reduces  to  a cross  sectional  area 
rather  than  volume.  The  effective  area  observed  by  the 
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mirror  at  D is  CB.  In  order  to  find  this  we  find  the 

equations  for  the  line  CB  and  AB,  and  intersect  them  to 

find  the  coordinates  of  C and  B. 

The  length  of  AB  = L 

AO  = OB  = L/2 


OD  = R 


Angle  DOB  = 0 


Equation  of  CB: 


+ L/2 


tanO 


tanQ 


Since  its  slope  is  -j./tan  because  CB  is  perpendicular  to 

OD.  Then,  finding  the  slope  for  AD,  we  get: 

RsinO 

m = 

RcosO  + L/2 


Equation  of  AD: 


RsinO 


— + -! 


(I./2)  RsinO 


RcosO  + L/2  RcosO  + L/2 

Intersecting  1 and  2,  we  get  the  coordimites  of  point  C: 


XC  = L/2 


(2Rsin  0 - R - (L/2)cos0) 
R + (L/2)cos0 


yc  = L 


cos  Rsin  0 - R - (L/2) cosO 
sin  R + (L/2)cos0 


Equation  of  the  line  BD  is 

RsinO 

y = } 

RcosO  - L/2 


L/2  sinO 
RcosO  - L/2 


Coordinates  of  point  B is  then: 
XB  - L/2 


yB  = ° 


J 


iiOHjitaBU’c 
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Using  the  distance  formula 


D 


we  get 


D = L 


cosG  Rcos0  + L/2 


sin0  R + (L/2)cos0 

At  this  point  because  of  the  relative  gain  factor  to  be 

used  we  can  take  the  problem  as  a volumetric  one,  then  the 

effective  volume  is 

V = nr2D 
e 

r being  the  radius  of  laser  beam. 

Relative  volume  ratio  is  taken  to  be  V_/V 

e'  max 

2 

where  V = nr  L 
max 

Therefore  the  volume  ratio  G(0)  is: 

^ /^\  _ cos  RcosG  + L/2 

b i \ U ) — 

sin  R + (L/2)cos0 


Assuming  R >>  L 


G1(0)  = 


2n 
cos  0 

sin© 


In  addition  to  this,  due  to  the  charging  of  the  parti- 
cle and  the  cross  sectional  effects  we  have  an  arbitrary 
correction  of  N(0)  which  is  the  number  of  light  particles 
(photons)  scattered  into  a unit  area  at  angle  0 [26] 

N (0)  = k/R2sin4 (0/2) 

K is  an  arbitrary  constant  which  depends  on  the  charge  and 
the  kinetic  energy  of  the  particle.  Therefore  our  correc- 


tion factor  becomes: 
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Fig.  A VI-1  The  log  intensity  versus  scattering  angle  graph  of  the  volumetric 
correction  factor,  treated  as  angular  gain. 


APPENDIX  VII 


SCATTERING  EXPECTATIONS 

There  are  a few  variables  that  enable  us  to  calculate 
the  expected  scattering  values,  and  compare  them  with  the 
real  data  points.  They  are  actually  the  basic  considera- 
tions of  our  experimental  system  design.  Aside  from 
several  other  factors  that  effect  the  measured  values  like 
the  coagulation  of  particles,  multiple  scattering  and 
variation  of  number  of  particles  in  the  scattering  region, 
an  approximate  estimation  of  the  first  order  can  be  done  by 
using  the  quantities  listed  below: 

1 - laser  output  power  and  bandwidth, 

2 - # particles  in  interaction  region  and  their  velo- 

cities , 

3 - The  angular  resolution  and  the  attenuation  of 

light  due  to  mirror  reflections  and  the  glan  air 
prism, 

4 - Photomultiplier  gain, 

5 - Digital  conversion. 

The  power  output  of  the  dye  laser  is  estimated  by  a 
comparitive  measurement  of  the  nitrogen  and  the  dye  laser 
outputs  (Figs.  8 a,b) . 

The  power  output  of  the  nitrogen  laser  is  calculated 
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by  the  readings  from  a thermopile  and  the  half  width  ob- 
tained from  the  oscilloscope  picture  through  a fast  photo 
diode. 

We  read  an  average  of  86  vi V output  from  the  thermopile 
at  10  pps,  counted  by  modulating  a frequency  counter. 
Calibrating  this  with  a He-Ne  laser  of  known  output  power 
which  gave  24  yV  at  2 mW  we  get: 


= 7.17  mW  average  output  power 

24 

7.17  m J/sec  1/10  pulse/sec  = .717  m J/pulse 
Half  width  = 6 ns 

. ^ _ .717  mJ/pulse  = 120  kW 

. . I - — g 

6 x 10  sec/pulse 


Comparing  the  amplitudes: 


P, 

_ dye  I 


'''dye 


Dye  laser  power  output: 

Tfl  U 

p = (fpTs)  '■*  X 120  = 81.98  kW 
15U 

Assuming  losses  of  45%  including  25%  due  to  the  Gian  air 
prism  that  is  used  to  polarize  the  incident  beam  and  20% 
due  to  the  focusing  lens 

1’  = 50  kW  approximately 

The  measured  half  width  for  dye  laser  output  is  also 
6 ns,  therefore  comparing  with  laser  we  have  approximate- 
ly .3  mJ/pulse 


Energy  of  a photon  = hC/A 


r 


Assuming  A = 627  nm 

-34  8 

F _ 6.62  x 10 joule  sec  x 3 x 10  m/sec  » 

1.6  x 10~19  3--u—  x 627  x 10~9  m 
sec 

Number  of  photons  per  pulse  is  then 

3 x 10  ^ joule/pulse 

1.6  x 10  '*'9  x 2 ev/photon 

14 

= 9.4  x 10  photon/pulse 

Assuming  1%  of  the  beam  is  scattered  equally  in  all  direc- 
tions and  knowing  the  solid  angle  covered  by  the  mirror: 


Solid  angle  = 


tt  (.63) 


7.4  x 10 


4 it  (36.6) 


the  number  of  photons  detected  per  pulse  is:  (taking  the 
number  of  photons  per  pulse  as  5 x 10^  due  to  losses) 

5 x 10  x8xl0~*-4x  10^'9  photons  detected 

per  pulse 

The  photomultiplier  specifications  are  given  in  table 
A VTI-1.  At  627  nm  the  radiant  cathode  sensitivity  E(A)  is 

E ( A ) = 30.4  mA/W 

10  -9 

with  4 x 10  photons  : we  get  12.8  x 10  joules 


. in10  ~ 
4 x 10  x 2 


= 12.8  x 10 


photon 


-19 

x 1.6  x 10  joule/ev  = 


Therefore  the  photocurrent  from  cathode  at  10  ns  will  be 


T , mA  12.8  x 10  W _ oo  n 

I . = 30.4  — x o 38.9  mA 

k W 10  x 10 


Taking  into  account  the  25%  attenuation  due  to  the 
analyser  prism  we  get  approximately  30  mA  output  from  the 


photomultipliers.  This  value  is  only  an  estimate  but  can 


be  used  to  determine  the  absolute  quantities. 


A better  approximation  can  be  calculated  by  knowing 
the  number  of  particles  in  the  interaction  region,  their 
sizes  and  shapes,  and  refractive  indexies. 

The  number  of  particles  in  the  interaction  region 
depend  or.  the  concentration  of  the  aerosols  in  their  sus- 
pension or  in  the  solvent,  the  liquid  flow  rate  and  the 
frequency  of  the  sinusoidal  voltage  at  which  they  are  gen- 
erated. These  are  the  variables  of  the  experiment.  They 
are  calculated  at  each  run  and  the  data  is  corrected  accord 
ingly . 

The  electronic  data  acquisition  system  lias  its  errors 
in  detecting  voltages.  These  are  given  in  Table  A VII-2. 
The  raw  data  is  corrected  due  to  the  relative  standard  de- 
viation of  these  functions. 


TABLE  A VII- 1 


Photomultiplier  Characteristics 


PMT  1 


PMT  2 


Type 

RCA7265 

RCA7265 

Serial  # 

P 57379 

P 12542 

Cathode 

Sensitivity 

220  yA/1 

255  yA/1 

Cathode  Red 

Sensitivity 

i-* 

o 

> 

1.25  yA 

Anode 

Sensitivi ty 

15500  A/1 

2700  A/1 

? 

2400  V 

2400  V 

Dark 

Current 

.045  yA 

.03  yA 

0 

193  OV 

220  V 

I 


APPENDIX  VIII 
DATA  TABULATION 

An  example  of  two  data  sets,  which  are  the  average  of 
5 consecutive  runs  each,  of  latex  particles,  are  given  in 
Fig.  A VIII-1  and  A VIII-2.  The  former  is  Channel  1,  and 
the  latter  is  channel  2.  The  incident  light  in  this  case 
is  polarized  perpendicular  to  the  plane  of  scattering.  As 
can  be  observed  from  both  of  these  graphs,  the  angle  at 
151.2  is  set  off,  and  the  axis  point  is  not  typed.  This 
data  point  is  considered  as  "out  of  range".  The  appear- 
ance of  the  out  of  range  point  is  because  of  the  miscoding 
of  the  paper  tape  in  one  of  the  five  runs  that  are  averaged. 
In  the  overall  averaging  this  point  is  discarded  as  well  as 
some  other  points  that  might  have  occurred  due  to  the 
scattering  from  coagulated  particles.  These  are  rejected 
if  they  are  within  the  range  of  the  error  bracket  set  by 

Mie  theory.  The  rejected  values  are  the  intensity  at  156.6° 

0 0 o 

of  Fig.  A VIII-1,  and  the  intensities  at  28.8  , 45.0  , 66.6, 
79.2°,  86.4°,  88.2°,  97.2°,  163.8°,  and  174.6  of  Fig. 

A VIII-2.  The  occurence  of  so  many  rejections  in  the 
second  channel  can  be  better  understood  by  observing  the 
graphs  of  the  two  channels  before  the  subtraction  of  the 
background  intensities,  which  are  given  in  Figs.  A VIII-3 
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and  A VIII-4. 

It.  is  clearly  seen  that,  while  the  first  channel  is 
giving  a smooth  output,  the  second  channel  is  very  incon- 
sistent. This  cannot  be  due  to  the  coagulation  of  the  la- 
tex particles  since,  the  scattered  light  is  detected  by 
both  channels  at  the  same  time,  from  the  same  orientation 
of  particles.  Therefore  we  conclude  that  it  is  a noise 
problem  of  the  second  channel  photomultiplier,  which  is 
random.  Nevertheless  the  data  points  corresponding  to 
this  effect  are  discarded. 

The  same  kind  of  reasoning  can  be  applied  to  the  clad 
osporium  and  NaCl  datas.  An  example  of  a set  of  cladospor 
ium  data  is  given  in  Figs.  A VIII-5  and  A VIII-6. 

Most  of  these  data  points  are  either  rejected  or 
omitted  within  reasonable  limits. 
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Fig.  A VIII-1  Channel  1 experimental  data  of  5 averag 
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Fig.  A VIII-2  Channel  2 experimental  data  of  5 averages. 

Incident  light  perpendicularly  polarized. 
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Fig.  A VIII-3  Data  of  Fig.  A VIII-1  before  the  subtraction  of 
background . 
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A VIII-4  Data  of  Fig.  A VIII-2  before  the  subtraction  of 
background . 


153 


«0  lOOJUftlfl’b 

«0 

'*iiWifO»*0 
10  lull*-*  to’i 

jo  nmm 

i)0  <uhs?vC6*I 
0»)  3uSw«»^g/#| 
X»*ld44?« 1 a*» 
0«  U'i»y6kfy*i 

SS  Mttubl 

fiO  1019  Oo«g#? 

oo  1i»ox??oyJ 

0J  ?u«»diOS.**4 

oo 

n unify- f 

o«  iciboots*? 
?l'  3»<>9»>,f.o*0 
00  lOSSI'O’j 

oo  lotouet* l 
oo  ^usoseolv 
00  30t»£»« if? 
oo  nor^ico*] 

00  lOOoBufSM 

oo  TuOe^ofo*! 

SS  J2WSSl:{ 

ss  ass*??  ir-\ 

00  OlUSOo’® 
oo  3w4*»9i  00  *£ 
00  SO/Jof®?*! 
OO  3v«>l<'0M  *1 
0 0 lo<>oC,iic,t  *1 
00  3yl£ft9bft*2 

oo  it  soao/ 1 
on  » I "s')  9 •? 

00  .1019(8/6*2 
20  H2b£Cyb.>-0 
Cn  Jo'.9S£''S,*2 

00  HoyiOi'f**? 
00 

00  Ho£„T2?9*2 


■>y‘,9S£''$,*2 
002192 Py* 


oot  * 


• *.0r 


20  T£')o-2M,"0 
JO  l,,iybyli*2 
CO  Aui  £ 

CO  ^ObyOiOi-? 

oO  SySiObPb’c 

00  /aOynt>"9*£ 

00  iy>i!f26i*£ 

88 

ss  wm \vi 

00  3v.'Se?£Po*£  1 

on  3unn£ini*i 

00  l.l2lOl)l£*{ 

oo  J02i92£ye 

00  lonSKSlfl 
00  3()/<'M£9*l 
00  3jS££o£i'£ 

00  Jj«t>b0n!,*£ 

on  ■>ilun,'(,{ 
oo  3U6 i22?o*0 
iiO  Wu?)."  2*t> 

00  Hyb2b.<BV'ft 
on  iCbyfsos-o 
00  3wn<infcC.8*6 
oo  .HO?0<U°2*ft 

oo  ho  l tUiC0*t> 

£ n HftC,u9iCl*0 
on  :H09entrf. 9 6 
00  njiyn^no'i 
00  WU?-<B>J*S. 
un  '^CiOtJV*, 
on  n .' w y n 2 n c • 9 
OP  30<’?Mi92*9 
yi  3j2IU(.6,9 
on  nv?r"£c }‘l 
£0  liJuoyi’O 

00  3Jb2f,0  * 

\ n H )"u9*'-o*T 

1 ’* 

2"  -49><sn\ *0 
20  nt,n\V£.»l>(, 

eo  3o/i>»B9t*o 

••••rtl-OO’O** 


-p 

0*9?l  0} 

e*»2i  *a 

»:e?t  -h 

2:S?  ° 

H ! c 

e:ij  M 

V * l I 

S:tt  • 

2-90  •"* 

? to 

lii  *-• 

g 

bbbbfSft  G 

S:tt  5 

S:SS  ■§ 

IW"*. 

»:2»  o 


e*J9  £.  r 

«A*t.s  n , 

8- SS  CXI 

8=11  S'i 

**os  ,« 

9*8*  S 

fi  * 9ft  (Of 

0*S*  <H 

2*£*  U 

9- 6? 

$WU  in 

Of  i 

9 *0f  £ 

18  § 

S-8  > 

6^65*12  < 

C • v.  1 

^’9! 

oOOQb*M  *y> 

cOftjr.Q*?* 

eooo**oi  ti, 

COOOJ'6  ^ 


154 


\ 


I0»10tuo9t>8*8 

lO-»xIMl46*8- 

o»  mi  i>,v?6*  i 
or  »r. tvi be x*£ 
|n*3x4it99n*9- 

10-ionuocaO’f* 
lO-mCl?.9i  *l« 
uc  mnnsco’S 
lo-9d?n<;i6»*9 
10--»j890V68*£ 

00  4.>»9ff  nT*T 

t0-10£(i989»*4 
oo  9x<iW?»'i 
in-iun:,?98 »*8 
00  1067/0(16*1 
in-3J09?e0D*9 
20-ms. 94?9*Z 

00  1x466094  *2 

ll  isssstn-l 

on  :iJVx?nll*t 
On  iu6Vi&B2*j 
o oust./ ?on* l 
o-3t£»92«i*» 

0-m49t691*9« 

u.iMMy-i- 

irMUim-r 

;o-icC*«2?*  *>**>- 

\0-3io6-?i/Ol*S- 
o*  nibbetfe 

ttzumivK 

on  tv/e oofo*2 

00  10;x£i?8*l 

KzmUWi: 

10-3x020919*9 
lo-iotbPoOtf  *n 

i8-3S5f?5r.?:l 

20  ?£li?£«/*0 
l o-l-onax i I /■ 
in-5Ji U£/ 6*2- 

1 o- mnsstfr  i • a- 

o'-  iooeiuo’J 
in  iuinb?i'(i 
to  •ix9r?T?n*i 
i n - 1 v o 0 n l / 0 f - 
\ d-luSt.  IlilM  *0- 

on  3j99iaf.fi 
l.n  3x£tV4T6*4 
on  1*2- 

00  3x69a2r«*l 
00  4xluU?0*4- 

oo  3xia?6CoM 
\o-3a2t>f9Cf  2 
00  msaOn/?*! 
on  309iorre*l- 

0 0 3x9 lim’2 
?0-Ju0/U°>’il* 
On  3tCi7 <°o*l 
00  1>ndO660*l 
?o-3c?c?£s;*e* 

I 0«9t??£CfO*I- 
in«3tS  n9nf>o*6 

OO  Sii»;on«I*2 
on  3,'Sf9i09*l 
on  3Sn9?5fi*» 

on 

«n 

00  , . 

in  v,M9Tfi*0 

on  t Oohtttu • ^ 

00  -x£  L\  i.*X  * l 

l"- -ill  < if*-!  *9- 
00  v*£ 

oo  io&o?ow*e 

s«  isfHsft-i 

on  ix948i(-9*8 
70  •</£!./!. <-4*0 
fln  1,.. nt  c;'f  i.  0 
..<•  !..x.6l'  t‘4 

On  3l>'9t>922*4 
00  1xl9/4«9  9 
on  :ivm94*'6!z 
in  :iu*,86»M>  9 
00  10064689*9 

»»»»Hl*00*0*h 


6*94 


.IfT,  3i  j .. 

aiiWiOPO*  1 

3ulwvttl( * t 


9 ’»Zt  ^ 
ooo5|*eii  a 


a 

w 

•H 


C-U 

R5 

9*99 

2*59 

• ’98 
9*99 

8*st 

2*18 

• *fc» 

9*4» 

8*8» 

0*»» 

*:sw 

oogoWct 

Vil 

|:« 

9*6? 

r-« 

?:?i 

9*5? 

9*91 

J*d 

8*f{, 

8:$$ 
ns 
r-ss. 

$:i*r 
M*£  S* 

9*16 

0*06 

e*w9 

66668*99 

r-\ s 

0*  »9 


*48 


319ST 


Fig.  A VIII-6  cladosporiun  data  of  channel  2 when  the  incident  light 
pendicularly  polarized. 
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